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ABSTRACT 



The analysis of vibration generation and transmission in shipboard 
machinery installations had defied rigorous mathematical description. 

After certain simplifying assumptions have been made, the vibration response 
of the installation can be predicted analytically. Theeobject of this thesis 
is to demonstrate the method by which the simplified system can be analyzed 
to give vibration trends and characteristics. The method incorporates the 
four pole parameter technique of analysis and a systematic study of the 
effects of installation modifications. 

The four pole parameters of a system are those expressions which 
inter-relate two basic input and output quantities. For vibrating mechanical 
systems these quantities are force and velocity. The four pole parameters 
must be derived or measured for the mechanical elements which transmit the 
force and velocity through the installation. In a shipboard installation 
these elements are the machine, mounting, foundation, hull, and water. 

After the four pole parameters are obtained for these five elements, 
matrix notation and combination permit a transmission study of the entire 
five element system. The computational assistance of a digital computer 
is an indispensable aid in this portion of the study. 

Numerous conclusions and recommendations were obtained from the 
analysis and from the results of the systematic study. The most important 
of these were: 

(1) For the purpose of vibration analysis a complex machinery 
installation can be represented by a relatively simplified 
system. This representation and the use of the four pole 
parameter technique are sufficient to establish trends and 
characteristics . 

(2) The stiffness of the resilient mounting has a dominant effect 
on the vibration transmission at all frequencies. 



ii 



(3) A Portsmouth mounting with decreased stiffness and increased 
damping and synthetic rubber damping applied to the foundation 
were the most effective installation modifications in providing 
overall transmission reduction. 

(4) As new ideas for reducing vibration transmission are intro- 
duced, they should be investigated by means of the four pole 
parameter technique. 

(5) A model of the approximate system of this thesis should be 
constructed to permit experimental verification of the four 
pole parameter technique. 



Thesis Supervisor: John R. Baylis 
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CHAPTER I 



INTRODUCTION 

1.1 Importance of Vibration Transmission 

The problem of noise control is one of major importance to naval 
engineers. A ship with a high noise level operates under the joint 
handicaps of being easily detectable by the enemy's listening devices and 
of reducing the effectiveness of its own sonar. Noise control implies 
the control of vibrations since acoustic or noise energy in a medium is 
synonymous with vibration transmission in that medium. A thorough under- 
standing of the vibration transmission characteristics of shipboard 
installations can provide the solution of the noise control problem. 

1.2 Current Developments 

Spurred by the Navy's need for quieter submarines, many facilities 
sind laboratories have been devoting time and effort towards finding a 
solution of the noise problem. Many studies are currently in progress 
seeking mecheinical means of reducing vibrations in machinery installations. 
Much of this is being financed by Navy contracts; thus, most of the major 
new developments have shipboard applications. 

The development receiving the most application at present in the 
Navy is the rubber-type resilient mount. These Navy mountings have been 
developed by the Engineering Experiment Station, Annapolis, Maryland and 
the Portsmouth Naval Shipyard, Portsmouth, New Hampshire ( 15 ) • The 
effectiveness of resilient mountings has been studied by Mr. A. 0. Sykes 
of the David Taylor Model Basin (13, 1*0. Studies have also been made of 
the properties of other rubber-like materials and of their effectiveness 
when used as mountings ( 17 )* 

It is becoming more apparent to interested engineers that the 
other elements in the vibration transmission path from source to water must 
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be investigated. Work has been done in determining vibration character- 
istics in machinery, the most common vibration source (12). Machinery 
foundation vibrations and techniques of damping them are being experimen- 
tally evaluated by the Sound and Vibration Section at the Electric Boat 
Company (20) . This work has emphasized the importance of the foundation 
in the vibration transmission path. 

The final item in the vibration transmission path is the ship's 
hull. Measurements of the vibration response of submarine hulls have 
been made. This empirical data gives a possible means of predicting the 
interacting effects of foundation, hull, and water (ll). 

It has been recognized that the time to apply noise control 
measures aboard ship is in the preliminary design and construction stages. 
This concept has been applied to the SSN 593> Thresher, to be commissioned 
this summer. A consistent set of specifications, setting allowable noise 
and vibration levels, have been developed for this ship. 

1.3 The System Under Study 

There are several means by which noise or vibration aboard ship 
can be transmitted from its source to the water. If an operating machine 
is considered as a typical noise source, its vibration can reach the hull 
and water by several paths . These are airborne transmission, supporting 
structure transmission, and transmission through connected piping. The 
vibration paths are best illustrated by the component block diagram in 



Figure I . 
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Figure I 

System Block Diagram 




The initial step in formulating an analytic approach to this 
transmission problem is to assume that each of these paths can be treated 
separately. This thesis is limited to an analytic study of the trans- 
mission through the supporting structure. The path under study, as 
indicated in Figure I, consists of machine, mounting, foundation, hull, 
and water . 

1.4 The Analytic Approach 

The exact describing equations of the transmission phenomena 
through the above mentioned mechanical elements would be an ideal analytic 
tool for the investigation of this problem. The actual system, however, 
defies any clear-cut description, and to obtain the describing equations 
is a Herculean task. The system is non-linear and possesses many degrees 
of freedom. The vibration of most machines is of a random nature. In 
addition, these machines are mounted on complex structures . 

To circumvent the above complexities, simplifying assumptions 
must be made to obtain a workable mathematical model. These assumptions 
are completely enumerated in Chapter II of this thesis. The simplified 



model is not intended to give the complete picture of vibration trans- 
mission in the actual installation. 

The object of this thesis is to make a systematic study of the 
analytically-described, simplified system. It is hoped that from this 
simplified study vibration trends and preferred design practice for the 
complex machinery system will become available . 
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CHAPTER II 
PROCEDURES 

2.1 The Preliminary Experiment 

The authors first became interested in an analytic approach to 
vibration transmission in shipboard installations while on temporary duty 
at the Admiralty Engineering Laboratory, West Drayton, England. A series 
of experiments were being conducted there to determine isolation charac- 
teristics of Admiralty resilient mountings. Results had been obtained 
from one of these basic experiments, and the authors investigated various 
approaches to the mathematical verification of this data. 

The experimental results studied were obtained by driving a beam 
supported by an Admiralty mounting with a typical vibration generator. 

This experiment is described in Appendix A-l. 

Most previous experiments investigating the characteristics of 
mountings treated the supported member as a lumped mass . This treatment 
facilitated use of a theory involving electrical analogies (lh). In this 
case the beam was used as the supported member because of its non -rigidity. 

A recent paper by R. D. Cavanaugh and J. E. Ruzicka provided a 
method by which a vibrating beam-spring system could be analyzed ( 4 ). 

This method was the four pole parameter technique . The four pole param- 
eter method had been widely used to analyze mechanical systems; therefore, 
it was selected as the technique to provide the mathematical verification 
of the beam-mounting experiment. 

2.2 The Four Pole Parameter Technique 
2.2.1 General 

There are many systems in nature which have two basic quantities 
as input and two basic quantities as output • An electric circuit element 
is one of these systems. In this case voltage and current are both input 
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and output. All these systems are called "four poles" since four quantities 
become inter-related by some action within the system. 

Electrical engineers have devised a technique of analysis for 
circuit elements called the four pole parameter technique. Each circuit 
element is considered independently, and complex circuitry is developed 
by connecting the various elements. Recently the four pole parameter 
technique has been applied to mechanical vibration problems where the 
mechanical element considered has a force and velocity as both input and 
output . 

A typical mechanical system can be represented in block diagram 
form as illustrated in Figure II . 

Figure II 
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F^ and are the input force and velocity, while F^ and are the 
output force and velocity. This system can be described in general by 
the following equations: 



F = e F + £ V 
1 11 2 12 2 



V 1 e 2l F 2 + 6 22 V 2 
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e . | 1 , € -| e 2 ]/ ari ^ e 22 are ^ our P°l e parameters of the mechanical 
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Using matrix notation, equations (l) and (2) become: 
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The determination of the four pole parameters comprises the major^ortion 
of the technique. Various rules must also be developed to permit combina- 
tion of the matrices to represent mechanical elements in various multiple 
configurations . 

2.2.2 Development of the Four Pole Parameters 

Basically the solution for the four pole parameters of the mechanical 
system can be reduced to the following procedure : 

1. Investigate the system and obtain its performance equations. 

2. Solve the equations imposing the proper boundary conditions. 

If non-linearities exist in the system it may be necessary to 
apply linearization techniques at this step of the procedure. 

3* Reduce the solutions to the form of equations (l) and (2). The 
four pole parameters can then be written by direct comparisons 
This procedure is applied to several basic mechanical elements in 
Appendix A-2. 

The above procedure can be specifically applied to the mechanical 
systems found in the analysis of vibration transmission. Equations (l) 
and (2) can be solved for the four pole parameters by imposing the condition 
that the output, which is station 2 in Figure II, is either blocked or 
free. The blocked condition implied the velocity at station 2 is zero; 
the free condition implies the force at station 2 is zero. 
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with station 2 blocked 
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The following definitions can be applied to equations (4) through 
( 7 ) to simplify the notation. 



(T™) = Force Transmission = •=— with station 2 blocked 

F F ± 



(T v ) - Velocity Transmission = rr— 
V V 1 



f V ? 
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Equations (4) through ( 7 ) become: 
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Equation ( 3 ) can now be written: 
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This development provides a convenient notation by which any 
system may be analyzed. 

2.2.3 Combination of Four Pole Matrices 

Multiple mechanical elements for which the four pole matrices 
have been derived can be connected in series. This is illustrated in 
Figure III . 



Figure III 
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The single four pole matrix [S ] is defined: 

s 

[S B ] = [S 1 3 X [Sg] X [S 3 J x [S n - 1] x [S n ] ( 9 ) 

Multiple elements can also be connected in parallel as illustrated 
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in Figure IV. 



Figure IV 




These n parallel subsystems can be represented by a single four pole 

matrix [S ]. 

P 




The derivation of the matrix [S^] is more complicated for the parallel 

than the series case. The solution for the four poles of [S ] is included 

P 

in Appendix A-3 • 

2.2.4 Application of Four Pole Parameter Technique to Preliminary Experiment 



From the preliminary experiment mentioned in section 2.1 and fully 
described in Appendix A-l, experimental results were available to test the 
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validity of the four pole parameter technique. Using this technique the 
isolation, which is the ratio of input force to output force for the system, 
was theoretically obtained as a function of driving frequency. The experi- 
mental results closely agreed with the values of isolation determined by the 
four pole parameter technique. This demonstrated the effectiveness of the 
technique * 

2.3 The General System and Analytic Study 
2-3-1 Description of the General System 

As mentioned in Chapter I, the overall system which is the subject 
of this thesis consists of the mechnaical elements in the structural path 
of vibration transmission from the machine source to the water. The 
block diagram of this system is Figure V. 

Figure V 




2*3-2 The Analytic Study 

The initial step in any analytic study is the enumeration of the 
simplifying assumptions. In this thesis these assumptions are: 

1. The transmission paths in Figure I can be considered separately. 
2- The principal contributors to structural vibration transmission 
are the force and velocity normal to the hull. All others are of 
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lesser importance . (This assumption has been corroborated by other 

investigators ( 17 ) • ) 

3 . The mechanical elements are elastic and linear. 

4. The vibrations in the system are periodic. 

Other assumptions have been made when analyzing the individual elements. 
These will be discussed in deriving the four poles of each mechanical 
element . 

It is recognized that the above assumptions will not permit an 
exact solution of the vibration transmission characteristics of an actual 
machinery installation. As stated in Chapter I, however, the purpose of 
this thesis is not to accurately predict vibration transmission for any 
particular mechanical system. Rather it is to provide a systematic 
study of a hypothetical installation such that trends caused by parameter 
variations can be predicted. A knowledge of these trends could indicate 
design procedures to reduce vibration transmission. 

With the use of the simplifying assumptions and the equations of 
dynamic analysis, the four pole parameter technique can be applied to 
the general system described in Section 2*3 .1. The dynamic analysis 
relates input force and velocity to output force and velocity for each 
element. The principal task is the determination of the four poles which 
appear in these equations. 

2.4 The Four Pole Parameters of the General System 
2.4.1 The Machine 

The problem of mathematically representing the vibration charac- 
teristics of an operating machine is a very difficult task. A study to 
evaluate the internal impedance of a typical machine has been made by the 
A. D. Little Co. (12). In this report an ideal machine is considered to 
possess many resonances and a complex response spectrum. Even though in 
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some cases the exciting internal forces are random, many machines rotate 
and thus exhibit a form of periodicity in their force generation. For 
this reason a harmonic force generator is considered as the vibration 
source in the A.' D. Little report. In their mathematical model a harmonic 
force generator drives a multi -mass-spring system in which the relative 
spring stiffnesses and mass sizes determine the resonances. Some investi- 
gators have concluded that this assumed lumped parameter behavior is valid 
at frequencies of only 150 cps or less (17). 

Recognizing the limitations of the lumped parameter treatment, 
this thesis postulates that the machine can be represented by a simple 
beam driven by a harmonic force generator. A simple beam by virtue of 
its structural properties possesses distributed parameters and eliminates 
the limitations imposed by the lumped parameter analysis. 

The four pole parameters for the beam are developed in Appendix 

A-4.1. 

2.4*2 The Mounting 

Of the mechanical systems in the structural transmission path, 
the mounting is the one best approximated by theory. Historically the 
installation of rubber-type resilient mountings was the first remedial 
step taken to reduce vibrations. Consequently, more mounting theories 
have been developed and tested than theories for the other mechanical 
elements . 

The theories used to develop the four pole parameters of the mounting 
in the general system are those proposed by D. V. Wright and A. C. Hagg 
(17, 7), A. 0. Sykes (13, 14), and J. C. Snowden (18). 

The first mounting considered in the general system is the standard 
Portsmouth mount, which is found in many shipboard installations today. 

The describing four poles for the mount are derived in Appendix A-4.2 
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according to the theory developed by Wright and Hagg. 

The second mounting considered in the general system is one 
constructed of a resilient rubber-type material similar to that bearing 
the trade name of Thiokol R. D. The four pole parameters for this mount 
are developed following the concepts of Sykes and Snowden in Appendix 
A-4.2 . 

2.4.3 The Foundation 

In most practical installations, machinery foundation structures 
are combinations of simple beams. This fact naturally suggests the use 
of beam vibration analysis to describe the foundation element. This 
analysis has been used by a number of investigators as a basic step in the 
description of the dynamics of foundation vibration ( 17 , 20 ). 

For the purpose of this investigation the foundation was simplified 
to be a beam resting on two columns. The describing four poles for this 
mechanical system are derived in Appendix A-4.3* 

2.4.4 The Hull 

The hull of a ship is a complex structure. Unlike other elements 
in the structural vibration transmission path, it cannot be represented by 
any single mechanical element or simple combination of elements. Experi- 
mental results are necessary to describe the response of typical hull 
sections ( 7 , 11 , 17 , 20 ). 

A mathematical model devised by D. V. Wright has provided a 
reasonable approximation of the vibration characteristics of a ship's hull. 
This model is used to obtain the four poles of the hull section in Appendix 

A-4.4. 

2.4.5 The Water 

The water is considered only in its effect as an attenuation 
element. Any entrained mass effect is included as part of the hull. The 
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water presents an effective specific resistance to the propagation of 
logitudinal pressure waves. This specific resistance is dependent on the 
density of sea water and the velocity of sound in that medium. 

The describing four poles are derived in Appendix A-4.5* 

2-5 Solution of the Problem 
2.5.I General 

The derivation of an expression for vibration transmission 
through the mechanical elements of the general system is the key to 
achieving the purpose of this thesis. The solution of this problem 
involves the combination of the five basic subsystems described in Section 
2.4 into a single complex system. The four pole parameters of the combined 
system are then obtained according to the equations of combination explained 
in Section 2.2*3* The vibration transmission is defined as a force ratio. 
The force transmitted into the water is considered the significant noise 
producing force in the system. Therefore^ the vibration transmission is 
defined as the ratio of the force into the water to the input force . 

Since many of the four poles in the system are functions of fre- 
quency^ the final ratio will be a function of frequency. Evaluation of ■ 
this ratio at sufficient frequencies will give the desired transmission 
versus frequency data. To conform to convention transmission is converted 
to decibels . 

The combination of the general system fourppoles and a sample 
calculation at one value of frequency are included in Appendix A-6. 

2.5-2 Computer Application 

This thesis could not have been completed without the computational 
assistance of the IBM 7^9 digital computer. The solution for vibration 
transmission at each frequency required a calculation of length comparable 
to that of Appendix A-6. This study required thousands of such calcula- 
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tions. The matrix notation of the four pole parameter technique was 
easily adapted to digital computer methods . 

A computer program was written to solve the combination of five 
four pole matrices. Provisions were made to change any of the mechanical 
elements in the system of Figure V. A complete sequence of this computer 
program logic is explained in Appendix A~7» A most convenient feature of 
the program was the use of the cathode ray tube oscilloscope to graphically 
present the results . The scope photograph eliminated the necessity of 
manual plotting. 

2.6 Systematic Study to Determine Trends 

With the selection of the general system of Figure V and the 
mechanical elements in Section 2.4 to represent this system, the problem to 
be considered was established. The four pole parameter technique described 
in Section 2.2 provided an analytical means to study the problem. The 
final phase of the study, however, was the most important. It was to 
devise a consistent, systematic method of changing the general system 
such that many different cases might be evaluated and compared. This 
comparison and evaluation formed the basis for determining the results, 
conclusions, and recommendations of this thesis. 

A detailed description of the systematic changes to the general 
system is given in Appendix A-5* 
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CHAPTER III 
RESULTS 

3.1 The Photo Plots 

The results of the analytic study of vibration transmission in 
shipboard installations are presented in the form of forty-one photo plots. 
These transmission versus frequency plots are included in Appendix B with 
one plot for each of the cases mentioned in Appendix A-5-7* 

3.2 Statements of Results 

The following statements of results are deduced from a study of 
the photo plots. To preserve the systematic procedure of the study, each 
case is examined separately. For detailed descriptions of the systems 
studied, refer to Appendix A~5 • 

Case 1 . This was the basic case involving System 1. It was 
selected as the basis for comparison since the elements comprising 
System 1 were typical of those presently found in shipboard 
installations . 

Case 2 . The System was the same as Case 1, except for the 
following : 

a. Mount damping increased. 

The first transmission peak was lowered. 

Case 3 . The System was the same as Case 1, except for the following 
a. Mount stiffness increased. 

The transmission over the frequency range studied was increased. 

Case 4 . The System was the same as Case 1, except for the following 
a. Mount stiffness decreased. 

The transmission over the frequency range studied was lowered by 
approximately 10 &b . 

Case 5 * The System was the same as Case 1, except for the 
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following : 

a. Thiokol R. D. mount replaced Portsmouth, mount. 

The first transmission peak was lowered while the second peak 
was removed. 

Case 6 . The System was the same as Case 1, except for the following: 
a. Hull damping was doubled. 

There was no change from Case 1. 



Case 7* 


The System was the same as Case 1 , except for the following: 


a. 


Mount damping increased. 


b. 


Mount stiffness decreased. 



The results were the same as Case 4 . 



Case 8. 


The System was the same as Case 1, except for the following: 


a . 


Thiokol R. Do mount replaced Portsmouth mount. 


b. 


Hull damping was doubled. 



The results were the same as Case 5- 



Case 9* 


The System was the same as Case 1, except for the following: 


a . 


Mount damping increased. 


b. 


Mount stiffness decreased. 


c . 


Hull damping increased twenty times. 



The transmission was two or three decibels below that of Case 7 at 
the lower frequencies . 

Case 10. The System was the same as Case 1, except for the 
following : 

a. Viscous foundation damping was added. 

The initial transmission slope was changed. The transmission was 
lowered at frequencies greater than 1700 cps. 

Case 11. The System was the same as Case. - 1, except for the 



following : 
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a. Thiokol R. D. and Portsmouth mounts in parallel. 

The transmission was increased at frequencies greater than 1000 cps. 
Case 12. The System was the same as Case 1, except for the 
following: 

a. Water resistance varied with frequency. 

The transmission was lowered at frequencies greater than 1700 cps. 
Case 13. The System was the same as Case 1, except for the 
following : 

a. Viscous foundation damping was added. 

b. Mount stiffness decreased. 

The first transmission peak was shifted from 21 cps to 11 cps. 
Transmission was decreased over the entire frequency range by 
at least 10 db. 

Case Ik. The System was the same as Case 1, except for the 
following : 

a. Thiokol R. D. foundation damping was added. 

Transmission was lowered appreciably at frequencies above 800 cps. 
Case 15. The System was the same as Case 1, except for the 
following : 

a. Thiokol R. D. and Portsmouth mounts in parallel. 

b. Portsmouth mount stiffness decreased. 

The results were the same as Case k to 500 cps. Above 500 cps, 
the results were the same as Case 11. 

Case 16. The System was the same as Case 1, e; cept for the 
following : 

a. Hull damping increased, maximum. 

The transmission was reduced over the entire frequency range. 

Case 17. The System was the same as Case 1, except for the 
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following : 

a. Thiokol R. D. mount replaced Portsmouth mount# 

b. Thiokol R. D# foundation damping was added. 

The first and second transmission peaks were removed. The 
transmission was reduced at frequencies above 800 cps. 

Case 18 # The System was the same as Case 1, except for the 
following : 

a. Thiokol R# D# foundation damping was added. 

b. Mount stiffness decreased. 

The vibration transmission was reduced 10 db# in the 10 to 800 cps 
range# Above 800 cps, the transmission was further decreased. 

Case 19 » The System was the same as Case 1 , except for the 
following : 

a. Thiokol R. D. foundation damping was added. 

b. Mount stiffness decreased. 

c. Mount damping increased. 

The results were the same as Case 18, except the first peak was 
removed. 

Case 20. This was the basic case involving System 2. Foundation 
mass increase was the only change from System 1. The transmission 
was lowered at frequencies from 10 to 1000 cps. The second and 
third transmission peaks were split . 

Case 21 . The System was the same as Case 20, except for the 
following : 

a. Mount stiffness decreased. 

The transmission over the frequency range studied was lowered by 
approximately 10 db. 

Case 22 . The System was the same as Case 20, except for the 
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following t 

a. Thiokol R. D. mount replaced Portsmouth mount. 

The first transmission peak was lowered; the second split peak was 
removed; and half of the third split peak was removed. 

Case 23 . The System was the same as Case 20 , except for the 
following : 

a. Mount damping increased. 

b. Mount stiffness decreased. 

c. Hull damping increased twenty tines. 

The first transmission peak was shifted below 10 cps, and 
transmission over the entire frequency range was lowered by 10 db. 
Case 24 . The System was the same as Case 20, except for the 
following j 

a. Viscous foundation damping was added. 

The split transmission peak at 3^-0 cps was lowered, and the 
transmission from 1000 to 10,000 cps was lowered by 5 to 10 db. 
Case 25 . The System was the same as Case 20, except for the 
following : 

a. Thiokol R. D. mount replaced Portsmouth mount. 

b. Thiokol R. D. foundation damping was added. 

The first transmission peak was removed while the second and third 
peaks were lowered. The transmission was considerably reduced at 
frequencies above 300 cps . 

Case 26 . The System was the same as Case 20, except for the 
following : 

a. Thiokol R. D. foundation damping was added. 

b. Mount stiffness decreased. 

The first transmission peak was shifted below 10 cps. With the 
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exception of the second peak, the vibration transmission was 
considerably reduced over the entire frequency range. 

Case 30* This was the basic case involving System 3* Foundation 
mass increase was the only change from System 2. The transmission 
was considerably lowered over the entire range except for peaks at 
800 and 4000 cps. 

Case 31. The System was the same as Case 30, except for the 
following : 

a. Mount stiffness decreased. 

The first transmission peak was removed, and the transmission 
reduced by 5 to 10 db. over the frequency range. 

Case 3^* The System was the same as Case 30, except for the 
following i 

a. Thiokol R. D. mount replaced Portsmouth mount. 

The first, second, and third transmission peaks were removed. 

Case 33 * The System was the same as Case 30, except for the 
following : 

a. Mount damping increased. 

b. Mount stiffness decreased. 

c. Hull damping increased twenty times. 

The results were the same as Case 31* 

Case ,3^ * The System was the same as Case 30, except for the 
f ollowing : 

a. Viscous foundation damping was added. 

The first transmission peak was split; the peaks at 800 and 2500 
cps were lowered. The peak at kOOO cps was removed. 

Case 35 * The System was the same as Case 30, except for the 



following : 
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a. Thiokol R. D. mount replaced Portsmouth, mount. 

b. Thiokol R. D. foundation damping was added. 

With the exception of the transmission peak at 200 cps, the 
vibration transmission was considerably lowered over the range of 
frequencies. 

Case 36 . The System was the same as Case 30, except for the 
following : 

a. Thiokol R. D. foundation damping was added. 

b. Mount stiffness decreased. 

The first transmission peak was removed, and the transmission 
reduced by 5 "to 10 db. up to 500 cps. Above 500 cps, the 
transmission was further reduced. 

Case 50 ♦ This was the basic case involving System k. The 
machine and foundation masses were ten times the values in 
System 1. This System was used as the basis of comparison for 
larger machinery installations. 

Case 51 « The System was the same as Case 50, except for the 
following : 

a. Mount stiffness decreased. 

The first transmission peak was shifted and lowered. The 
transmission was reduced by 5 to 10 db. over the 
Case 52 . The System was the same as Case 50, except for the 
following t 

a. Thiokol R. D. mount replaced Portsmouth mount. 

The first transmission peak was removed, and the transmission was 
lowered by 5 to 10 db. over the entire range. 

Case 53 - The System was the same as Case 50, except for the 



following : 
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a. Mount damping increased. 

b. Mount stiffness decreased. 

c. Hull damping increased twenty times. 

The first transmission peak was considerably lowered. The 
transmission was lowered over the entire range by approximately 
10 db. 

Case $4. The System was the same as Case 50 , except for the 
following : 

a. Viscous foundation damping was added. 

The first transmission peak was lowered; the third peak was 
eliminated. The transmission was lowered by 5 to 10 db. from 1000 
to 10,000 cps. 

Case 55* The System was the same as Case 50, except for the 
following : 

a. Thiokol R. D. mount replaced Portsmouth mount. 

b. Thiokol R. D. foundation damping was added. 

The first transmission peak was removed, and the transmission was 
reduced by 10 to 15 db. to 300 cps. Above 300 cps, it was reduced 
by 20 to 40 db. 

Case 56 . The System was the same as Case $ 0 , except for the 
following : 

a. Thiokol R. D. foundation damping was added. 

b. Mount stiffness decreased. 

The first transmission peak was lowered and shifted to 12 cps. 

The second peak was reduced by 35 db. At higher frequencies the 
transmission was extremely low. 

Case 57* The System was the same as Case 5 0 > except for the 



following: 



(25) 



a. Thiokol R. D. foundation damping was added. 

b. Mount stiffness decreased. 

c. Mount damping increased. 

The first transmission peak was removed. The second peak was 
reduced by 35 <3b . At higher frequencies the transmission was 
extremely low. 
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CHAPTER IV 

DISCUSSIOH OF RESULTS 

4.1 Preliminary Discussion 

The most important single result of this thesis is the technique: 
of analytic description of a complex shipboard machinery installation. 

This technique formulates a workable analytic tool for the study of 
vibration transmission characteristics of these complex systems. In 
essence this technique is the four pole parameter method. It is recom- 
mended that this method of - analysis be considered in future design inves- 
tigations . 

4.2 Systematic Changes 
4.2.1 Mounting 

The properties of the resilient mounting had a profound effect 
on the vibration transmission. A study of the photo plots showed that 
stiffness in the Portsmouth mounting had a dominant iffect in the cases 
in which it was altered. When stiffness was increased, the transmission 
at all frequencies was increased; when stiffness was decreased, the 
transmission at all frequencies was lowered by 10 db. It is recommended 
that the design engineer make stiffness as low as possible and yet 
consistent with shock, alignment, and other factors influencing mount 
design. 

Damping in the Portsmouth mounting was not nearly as influential 
as mount stiffness. However, increased damping was Effective in decreasing 
the transmission peaks at lower frequencies. It is recommended that an 
effort be made to increase! the damping of the Portsmouth mount by increas- 
ing the material damping or by introducing an external damping device in 
parallel with the mount. 



The replacement of the Portsmouth mounting by a synthetic rubber 
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mounting of Thiokol R. D. material showed interesting results. The 
Thiokol R. D. mount lowered the first transmission peak and eliminated the 
second. Therefore, the use of the Thiokol R. D. mount appears desirable 
in systems with high noise transmission at low frequencies or in systems 
in which low frequency noise suppression is paramount. The stiffness of 
Thiokol R. D. increases exponentially with frequency (l8). This decreases 
the effectiveness of Thiokol R. D. in suppressing vibration transmission 
at higher frequencies. The development of a material to correct this 
undesirable trend of Thiokol R. D. is recommended. 

4 .2 .2 Foundation 

Of the various foundation parameters the mass most affected the 
vibration transmission. Three different foundation masses were investi- 
gated with the result that vibration transmission was lowered as the 
foundation mass increased. Within the limitations of weight and space, 
foundations of vibrating machinery should be made as massive as possible. 

Viscous foundation damping and damping utilizing the Thiokol R. D. 
material were studied. Foundation damping was effective at frequencies 
above 1000 cps for the viscous type and above 800 cps for the Thiokol R. D. 
type. However, the vibration transmission was reduced more by the Thiokol 
R. D. type. 

Foundation damping is recommended in installations where high 
frequency noise transmission is a problem. 

4.2.3 Hull 

The effects of increasing the hull damping first became apparent 
in the mathematical model used when the damping was increased by an order 
of twenty. As hull damping was further increased, the vibration transmission 
was reduced at all frequencies. 

These results are at best qualitative. In this model lumped 
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parameter damping was used, and no attempt was made to study the effects 
of distributed damping. In practice the latter would be the type applied 
to the hull. A study to quantitatively evaluate the effects of distributed 
hull damping would be useful in the extension of an analytic study. 

4.2.4 Combined Changes 

In a series of cases studied the Portsmouth mount Stiffness was 
decreased, the mount damping increased, and b[i1 1 damping increased twenty 
times. This was done with the four systems studied. In general, as the 
masses in the systems were increased, the vibration transmission was 
reduced. The results of this series study were compared with those cases 
in which mount stiffness alone was changed. This comparison showed that 
the change; in mount stiffness overrode the changes in mount damping and 
hull damping. 

A single case of decreased Portsmouth mount stiffness and increased 
mount damping showed that the vibration transmission was lowered by exactly 
the same amount as the case in which stiffness alone was decreased. This 
further supported the conclusion that mount stiffness was one of the most 
important factors to be considered in any vibration transmission study . 

When a mount of Thiokol R. D. was placed in parallel with a 
Portsmouth mount, the transmission reduction properties were not improved. 
The addition of the Thiokol R. D. mount actually increased the transmission 
at frequencies greater than 1000 cps. This case was again studied with 
the Portsmouth mount stiffness decreased. Results showed that this case 
was not as good as that in which the single Portsmouth mount with decreased 
stiffness was examined. Therefore, it is concluded that placing the 
Thiokol R. D. mounting in parallel with a Portsmouth mounting would not be 
advantageous in reducing vibration transmission. If parallel mountings 
are to be used, as recommended by Snowden in Reference (l8), care must be 
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used in the selection of the high damping synthetic rubber material. 

In another series of cases the Portsmouth mount was replaced by 
a Thiokol R. D. mount and Thiokol R. D. damping was applied to the founda- 
tion. This was done with the four systems studied. As the masses in the 
system were increased, the vibration transmission decreased markedly. In 
any one system the combination of Thiokol R. D. mount and Thiokol R. D. 
foundation damping was most effective in reducing transmission. It is 
recommended that this combination be investigated in an experimental 
installation. 

A Portsmouth mount with decreased stiffness and Thiokol R. D. 
foundation damping comprised the changes examined in another series of 
cases. Again this was done with the four systems. The vibration 
transmission was significantly decreased as the masses in the system were 
increased. Within a particular system this combination was less effective 
at frequencies of less than AO cps than the Thiokol R. D. mount and Thiokol 
R. D. foundation damping combination. However, it was superior at 
frequencies higher than 40 cps. This combination should also be investi- 
gated in an experimental installation. 

A Portsmouth mount with decreased stiffness and increased damping 
and Thiokol R. D. foundation damping comprised the changes in the final 
series examined. The only transmission change^ from the previous series 
was the elimination of the first transmission peak. 

^•3 General Considerations 

The results of the systematic study largely verified existing 
concepts in reducing vibration transmission in shipboard machinery 
installations. This would indicate that the assumptions and the four pole 
parameter method of analysis were sufficient to establish the vibration 
trends and characteristics proposed at the beginning of this thesis. 
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As new ideas for reducing vibration transmission are introduced, it is 
proposed that they be investigated by means of the technique of this 
thesis. 

The four pole parameter technique enables use of empirical data 
obtained independently for each mechanical element in the vibration path. 
Empirical data can be used for elements which defy analytic description. 
It would appear that the use of empirical data for the element four poles 
would more closely approximate the vibration characteristics of an actual 
system. 

An interesting extension of this thesis would be an optimization 
study to determine the best combination of elements in the vibration 
transmission paths. Since shock specifications also govern the element 
selection, the optimization must include shpck as well as vibration 
reduction considerations. 

The construction of a model of the approximate system, as 
represented in Figure A~3 of Appendix A, would permit an experimental 
verification c£ the four pole parameter technique. 
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CHAPTER V 
CONCLUSIONS 

As logically developed in Chapter IV the following statements 
are listed as the conclusions of this thesis: 

1. The four pole parameter technique provides a sufficiently 
powerful method of analyzing vibration transmission. Furthermore, 
the four pole parameter technique is easily adapted to a digital 
computer program. 

2. For the purpose of vibration analysis a complex machinery 
installation can be represented by a relatively simplified system. 
This representation is: sufficient to establish vibration trends 
and characteristics. 

3. The stiffness of the resilient mounting has a dominant effect 
on the vibration transmission at all frequencies. 

k. The only effect of increased damping in the Portsmouth mount 
was to lower the first vibration transmission peak. 

5 . A mount fabricated from a high damping synthetic rubber such 
as Thiokol R. D. is highly effective in reducing vibration 
transmission at low frequencies. The same mount is not as effective 
as the Portsmouth mount at higher frequencies. 

6. As the mass of the foundation is increased, the vibration 
transmission is lowered. 

7. Foundation damping was effective at frequencies above 1000 
cps for the viscous type and 800 cps for the Thiokol R. D. type. 

The Thiokol R. D. type was more effective. 

8. Increasing hull damping decreases vibration transmission at '.all 
frequencies . 

9 . Placing a Thiokol R. D. type mounting dn parallel with a 
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Portsmouth mount as suggested in Reference (l8) actually had a 
detrimental effect on vibration reduction. 

10. Replacing a normal mount by a Thiokol R. D. mount and applying 
Thiokol R. D. damping to the foundation were most effective in 
reducing vibration transmission. 

11. Decreasing the stiffness of a Portsmouth mount and applying 
Thiokol R. D. damping to the foundation were superior to the 
system mentioned in the previous conclusion in reducing vibrations 
at frequencies above 40 cps. 

12. A Portsmouth mount with decreased stiffness and increased 
damping and Thiokol R. D. damping applied to the foundation were 
the most effective combination in overall transmission reduction. 
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CHAPTER VI 
RECOMMENDATIONS 

As stated or inferred in Chapter IV the following are the 
recommendations of this thesis: 

1. The four pole parameter method of analysis should be considered 
in future design investigations involving vibration transmission. 

2. As new ideas for reducing vibration transmission are introduced, 
they should be investigated by means of the technique of this thesis. 

3. Empirical data should be used for element four pole parameters 
to more closely approximate the vibration characteristics of an 
actual system. 

b, An optimization study to determine the best combination of 
elements in the vibration transmission paths should be conducted. 

It should include shock as well as vibration considerations. 

A model of the approximate system of this thesis should be 
constructed to permit experimental verification of the four pole 
parameter technique . 

6. Mount stiffness should be as low as possible and yet consis- 
tent with shock, alignment, and other factors influencing the 
design. 

7 . The material damping of the Portsmouth mount should be 
increased, or an external damping device placed in parallel with 
the mount. 

8. The Thiokol R. D. mount should be used in systems with high 
noise transmission at low frequencies or in systems in which low 
frequency noise suppression is paramount. 

9* A high damping mounting material should be developed in 
which stiffness remains constant or decreases with frequency. 
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10. Within the limitations of weight and space, foundations of 
vibrating machinery should be as massive as possible. 

11. Foundation damping should be utilized in installations where 
high frequency noise transmission is a problem. 

12. A study to quantitatively evaluate the effects of distributed 
hull damping should be conducted. 

13* Caution should be exercised in' the selection of a high damping 
synthetic rubber mount to be placed in parallel with a standard 
mount. 

l4. The effects of replacing a normal mount by a Thiokol R. D. 
mount and applying Thiokol R. D. damping to the foundation should 
be investigated in an experimental installation. 

15* The effects of decreasing the stiffness of a Portsmouth mount 
and applying Thiokol R. D. damping to the foundation should be 
investigated in an experimental installation. 

16. The effects of decreasing the stiffness and increasing the 
damping of a Portsmouth mount and applying Thiokol R. D. damping 
to the foundation should be investigated in an experimental 



installation . 
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APPENDIX A 



DETAILS OF PROCEDURE 



A-l The Beam-Mount Ing System 
A-l . 1 Exper iment 
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This experiment was conducted to determine the isolation charac- 
teristics of a typical Admiralty mounting. The experimental system consisted 
of a fifty pound beam driven by a harmonic vibration generator. The beam 
was supported by the mounting which rested on a "rigid" foundation. The 
components and instrumentation are illustrated in Figure A-l. 

The foundation system was not, in fact, rigid. However, its natural 
frequency was much lower than the rigid body mode frequency of the beam- 
mounting system. This insured no resonance interaction, and the foundation 
mass, nip, in Figure A-l, responded in a rigid manner. The vibration 
generator imparted a velocity to the mass nip. The velocity into the 
mounting of compliance and the velocity at the free side of the beam 
were picked up by the barium titanate accelerometers, A^ and A^ respec- 
tively. The isolation in this case is the ratio of to V^. Both this 
and the driving frequency were recorded, and the isolation was converted 
to decibels. This was repeated over the frequency range of 10 to 3000 cps. 

To obtain a more realistic model in which the beam was driven 
from above, as a machine normally would be, use of the transmissibility 
theorem was necessary. This stated that the transmission (or isolation) 
in one direction for linear elastic systems is equivalent to the force 
transmission in the opposite direction: 




This new system could be represented: 
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A,,A t - Barium T»tanate 
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rr» t ~ Mass of B«aM 

m f - Mass of Fouhoatioih 
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The physical constants of the system of the experiment were as follows 
Beam: Length: L = 60 in. 

Thickness: t = -75 in. 

Weight: = 50 lb. 

.D 

Frequencies : 




g 

where : E - 30 x 10 psi 



W t 

I - moment of inertia. = -r— - — — 

12 L 



^ = mass density 
Describing Functions: 



JZL (^f) and ^ given in Reference (l6) 
Mounting : Compliance : ^ 

A-1.2 Theory 

The beam-mounting system was represented schematically: 
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[S ] was the matrix of four pole parameters for the beam and 
[ ] for the mounting . 

Then: 




or by individual steps : 




The matrices in terms of their four poles were expressed: 







1 


i 






H 

H 

w 


£ 12 


S 1 




i 


1 






1 

ro 

H 


e 22 



( 1 ) 



(la) 



(lb) 



( 3 ) 



\\\\w 



and 



(bo) 



s ? 




2 

e ll 


2 

£ 12 


°2 




2 


2 






1 

ro 

H 


<T> 

ro 

ro 

i — 



w 



As developed in Section 2.2.2, the general four pole parameter 
matrix is: 
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In this proof of the experiment the mounting was treated as a 
simple spring. From Appendix A~2, the four poles for a simple spring are: 
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The compliance of the spring is C^, while the mobility is 

The beam was analyzed following the techniques of Ruzicka and 
Cavanaugh ( 4 ). With negligible compression assumed in the transverse 
direction and from the transmissibility theory, the four poles of the beam 
were evaluated : 

1 b 1 f 
^tT^ = ^tT) = 1 



(m 21 ) d - o 



(i — ■) = iL— 
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where was the mobility at the center of the beam. This was taken 

from Reference ( 4 ): 
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The matrix [ ] for the beam was then written : 
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Equations (5) and. ( 7 ) were substituted, in equation (l): 
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To check the experimental results a plot of isolation I 9 which is 
F i 

defined as — > versus frequency was desired. Equation (8) was rewritten: 

F 3 



V. 



jwc 

7l 

“ 2 ' “'2 

\ = (j^) f 3 + v 3 



F i = (l + f~) F + -2=- 

1 M(^) 3 M(£) 



Station 3 in "the experiment was a rigid, support; therefore, V 3 = 0. 
The isolation was obtained.: 






M(|) 



The mobility of the beam is expressed in Equation (6). 

A-1.3 Comparison of Experiment and Theory 

The curves of isolation in Figure A-2 show that the four pole 



parameter technique provided excellent theoretical verification of the 
experimentally measured values of isolation. These results indicated the 
value of the four pole analysis and suggested the use of it in the study 
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FIGURE A-2 

PRELIMINARY EXPERIMENT 




(QP) N0I1V10S I 



FREQUENCY (cps) 



( 43 ) 



of more complex systems. 

A-2 Development of Four Poles 
A-2.1 Mass Four Poles 




The performance equations are 



V = V 
1 2 



because the mass is a rigid body and 



dV. dV 

F 1 " F 2 = m dt“ = m dT 



( 1 ) 



( 2 ) 



from the application of Newton's Law. 

Since the applied force is usually considered as harmonic, the 
other force and the velocities will also be harmonic functions . 
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and V 2 q are time 



Substituting the expression for Equation ( 2 ) becomes: 
F x = F 2 + (>m) V 2 

while (l) remains 



V = V 
1 2 



By inspection the four pole parameters are: 



£ 11 = 1 ' e l2 = € 21 = °' 311(1 £ 22 = 1 * 

A-2 .2 Spring Four Poles 



k 




The performance equations are 



F = F 
1 2 



( 1 ) 



because the spring possesses no inertia and 

F^ = k [x^ - Xg ] where x is displacement. (2) 

The velocity is the time derivative of the displacement . 

From the harmonic velocity expressions of the previous section: 
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Equations (l) and (2) are rewritten: 
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By inspection the four pole parameters are ; 
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A-2.3 Damper Four Poles 
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The performance equations are 




because the damper possesses no inertia and 

F 1 = r(V l " V 

from the definition of damping. 

These are rewritten 



F = F 
1 2 



V = — F + V 
1 r 2 2' 



By inspection the four poles are 
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A— 3 Parallel Matrix Combination 



( 1 ) 

( 2 ) 



If n systems are connected in parallel then a single four pole 



matrix can be written connecting the common input and output points. If 



this matrix is written [S ], then 
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Let the r matrix, [S^], be given by: 
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where r varies from 1 to n. Then, in (l) above 
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A more complete discussion is given in reference (6). 

A-4 Development of Four Poles of General System 
A-4.1 The Machine 

As stated in Section 2.4.1 the machine was represented by a 
uniform elastic beam driven by a harmonic force generator. This system 
was selected because it was multi-resonant and possessed distributed 
parameters . 

The four poles of this system are those briefly described in 
Appendix A-1.2 for the beam in the preliminary experiment. These four 
poles were = 1, ~ = 0, and = 1 . From the method 
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of Reference (4) , the beam mobility M(— ) will be derived. 



The basic equation of free vibration is; 



tl 2-%- A £-2- = 0 

d x d t 

where: I = moment of inertia of beam's cross section 

E = modulus of elasticity 
= mass density 
A = beam cross-section area 
y = transverse displacement 
x = distance along length of beam 
L = length of beam 
t — time 

Assuming a harmonic solution: 

y = fMe** 

If this is substituted, equation (l) becomes: 
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and W = frequency of vibration. 

The general solution of (3) is: 

f(x) = (sinh^x + sin^x) + (cosh^x + cos^x) + ( sinh 

^x - sin ^x) + (coshy5x -tcos^x) (4) 



T^e boundary conditions for a free-free beam are: 

a g f( o) Q 

dx 2 



( 5 ) 



(») 



d g f(L) 

dx 2 

,3 



= o 



a 3 f(o) _ 0 



dx" 



d 3 f(L) _ 



= 0 



dx*“ 



( 6 ) 



(7) 



( 8 ) 



Equations (5) and (6) require no shear forces at the ends of the beam; 
equations (j) and (8) state that there are no end moments. 

When the boundary conditions are applied to equation (4), the 
characteristic frequency equation is obtained. 

cos /^n ^ cos h /$ n L = 1 (9) 

This is a transcendental equation. Corresponding to each eigenvalue 

will be a value of f^(x) . The term f (x) represents the mode shape in 

which the beam vibrates at the frequency defined by^^. The relationship 

between f (x) and A is: 
n' ' rn 

fi^x) = C 2 [cosh^S n x + cos^ - o< n ( sMh^ n x + sin^ x) ] (10) 



where 



sinh A L + sin A L 

^ = ,rn / n 

to. cosh L - cos A L 



r n 



( 11 ) 



C 2 = arbitrary constant 

Define : • ) 

f (x) = C d) (x) . 
n x 2'n 

Then in equation (10): 

(J? n (x) = cosh x + cos x (sinh x + sin x) 

The term (J> n (x) will satisfy the equation: 



(^9) 



dty) 
' 1 



dx 



n „ 4 i 



r L 

L4A dx = 0 n 



It is a normal function which forms an orthogonal set in the interval 
O^xiL. The condition of orthogonality: 

► L 
'0 

However, when n = m: 

J^>n (x) ( L )i 2 

D. Young and R. Felgar have evaluated these (j) functions for beams with a 

variety of boundary conditions (l 6). 

With the values of (J> n and the corresponding values of UJ^, the 

forced vibration motion is considered. 

iwt 

If a harmonic force F^e 0 is applied to the beam, the displacement 
according to reference (l) is: 



y(x) = 



F o e 



jwt 



go 



“b 



<*> 

Z. ^2 2 + b°dy motion 



( 12 ) 



n = lu> -u) 
n 



The rigid body motion is f (x) e^*^. The velocity is the derivative of 
the displacement: 

. ju/F e Jwt <T“ 2 ( x ) 

y(x) = Z “§ J + jwf Q (x) e 0 * 



“B 



n = 1 (J ~u) 
n 



( 13 ) 
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From the definition of mobility at the driving point p: 

v(p) _ y(p) 



M = 



F(P) F 
0 



c?o 



M(p) 



= 2_ 



4n (p) ^0 



"b ^t-iw 2 - w 2 f o 

n 



( 1*0 



As a rigid body the beam acts like a mass and possesses the mobility of 
this lumped parameter. From Appendix A-2.1: 

Mass Mobility - — ^ — — 

Therefore, in equation (14): 
jwf . 



0__lL 



0 



V * 



The final expression for beam mobility is: 



CSO 



M (p) y~ 

“b fr-i u 2 -cj 2 “b 1 

n 



(15) 



Since the beam in the machine representation was considered to be 
driven at its center, the four pole e ^ required the inverse mobility 
at the center of the beam. This was expressed from equation (15) as: 

’ *h.W-u? V*' 



M(: 



n 



The machine could now be represented in terms of its four poles. 

A— 4 .2 The Mounting 

According to Section 2.4.2, two different mountings were studied. 
The first was the standard Portsmouth mounting which was analyzed following 
the method of Wright and Hagg (17)> 

The empirical results of tests on actual mountings were used by 
Wright and Hagg. The mount characteristics were approximated by use of the 



( 51 ) 

electrical impedance analogy, one used extensively in vibration analysis 
(8). The mounting was represented by a shunt impedance consisting of real 
and imaginary parts. In the language of the analogy, impedance is defined 
as a force to velocity ratio. 

*=£ ( 1 ) 

In this case input impedance was considered for the approximation. 
From the analogy mechanical damping corresponds to electrical resistance 
or the real part of the impedance. Mass and spring-like properties 
correspond to electrical reactance or the imaginary part of the impedance. 
Thus the approximation included all the important lumped parameters . 

These parameters, damping constant R and spring stiffness K, were 
established to be : 



R _ K 

" 2*fQ Q 

k « 2 W£?(l + f-) 

- k — “ 

g(l + T- 2 -) 
t 

where : 



( 2 ) 

( 3 ) 



fg = rated natural frequency of the mount at its rated load, cps, 
empirically determined. 

Qq = maximum transmissibility at f^, dimensionless, empirically 
determined. 

f^ j= frequency of the first high frequency transmission peak, cps, 

empirically determined. 

W 

— = rated mass load on the mount, lbm, empirically determined. 

S 

f = driving or excitijag frequency, cps. 

The shunt impedance Z is therefore given by: 
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Z = R + 4- (4) 

3 “ v ' 

The impedance analogy indicates the following equations hold for 

a shunt impedance: 




F = F 
in out 

F 


(5) 


„ out . . , 

Z - ■, open circuit 

V in 


(6) 


V. = V , , short circuit 
m out' 


(7) 


Equations (5) to (j) are translated 
notation as in Section 2 .2 : 


into fQiir pole parameter 


F = F 
1 2 

„ F 


(8) 


1 - z -2. 
"W” 


(9) 


station 2 free 


(10) 



Therefore, the four pole parameters for this standard Portsmouth 
mounting are 

€ 11 = ^ € 12 = °' e 21 = b ^ e 22 = X * 

Combining equations (2), (3)j and- (4) the following expression for 



is obtained. 



(53) 



= wQ 0 
€ 21 K 



1 + jQ 0 

, Q o + h 



( 11 ) 



where co is the forcing frequency, rad / sec . 

The second mounting considered was treated similarly except that 
different expressions were used for R and K. These were derived by A. 0. 
Sykes (l4): 



. 2rG 

uJ(l - r 2 ) 




( 12 ) 



K - 0(f) 



(13) 



where : 



S = mounting cross-sectional area 
L = mounting length 
G = mounting elastic modulus 
r = mounting damping factor 
b-' =■ forcing frequency 

The elastic modulus and dancing factor were given as empirical 
results for a synthetic high damping rubber called Thiokol R. D. by J. C. 
Snowden ( 18 ) . 

Approximate formulae were determined for G and r from the empirical 



data: 



G = 72.6 OJ 



. 0-4 



0.325 



(lb) 

( 15 ) 



r = 0.1 U 

Thus R and K could be represented as functions of frequency and the new four 



pole parameter e ^ evaluated. 



W 



' 21 0(f) (A 2 + 1) 



[ A £ j ] 



( 5*0 



where 




A-4.3 The Foundation 

As stated in Section 2.4.3, the foundation was considered a beam 
resting on two columnar structures . The two columns presented two possible, 
identical paths for vibration transmission. Since the paths were identical, 
only one was studied explicitly. 

The four poles of the foundation structure were a series combination 
of the individual four poles of a beam and column. 



e n e i 2 




1 m(l7 




1 0 


e„, e 




0 1 




M_ 1 


|_21 22j 








L c J 



The combined element four poles are 

M C 1 

6 11 “ 1 + M^Lp e l2 ~ M(Lp e 21 “ M c' 3X111 e 22 = 1 * 

The term M(L) is the transfer mobility of the beam relating the 
driving point with the junction of the beam and column structure. Initially, 
this beam was considered as end-supported for the purposes of deriving its 
mobility. Correspondence with the Admiralty Engineering Laboratory, where 
an experiment driving an end-supported beam was being conducted, indicated 
that for vibration analysis the end-supported beam behaved more like a 
free-free beam. For this reason the foundation beam was treated much as 
the beam in Appendix A— 4 . 1 . The expression for beam mobility is: 

”b fr-1 w 2 -w 2 v* 

n 

where the (j) ' s are the normal functions at the driven point (^-) and the beam 
column junction (L) . Again these functions were taken from Reference (l6). 
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From Reference (l) the expression for the lateral deflection of 
a column being driven by a harmonic force F^e is: 



n-1 



y(x) » f e jut + rj 

n = 1 r 




■where : 



mass density 
L - column length 



A = column cross; section area 



The column mobility can be derived using the same methods used for the 



beam in Appendix A-4.1. The velocity is: 




n-1 



n = 1 



2(-l) 2 

(«£ -“ 2 > 



Now the definitions of mobility and the rigid mass mobility are applied. 
From this the mobility of the column is obtained: 



A-4 . 4 The Hull 

As stated in Section 2.4.4, the work of D. V. Wright was consulted 
again for a workable analytic representation of the ship's hull (7)* 

Wright considered the hull a series impedance, which he represented 
as a result of analyzing empirical data obtained from typical hull impedance, 
measurements (11, 20) . As in the case of the mounting the impedance 
consisted of a real and imaginary part, the resistive and reactive properties. 
These were represented as complicated quantities dependent on nine separate 
lumped parameters . 




n-1 



n 
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R = R 



s 



X = X 



s 



R 2 

+ R x f + y ~ + 



# f K 

TCWS S 

+ 193 " 2 rtf 



R 

-P. 



1 + R ( 



,2 ,193 2*f 



rtW f 
P 



K 



) 



R 



/ 193 _ 2rtf\ 

\ -p k ' 



kW f 
P 



1 + R ( 



,2 ,193 2rrf 



rtW f 
P 



K 



) 



( 1 ) 



( 2 ) 



Z = R + jX (3) 

where R , R n , R^, R , K , K , W , W are the nine lumped, parameters and 
s' 1' 2’ p' s' p' s' p 

f is the driving frequency. 

The nine parameters were considered circuit elements in the 
impedance analogy, making up the series impedance Z. 




The impedance analogy gives : 



z = 


in 

V. 


(4) 




m 




V in 


= V 

out 


(5) 


F. 

in 


= F ,, open circuit 
out 


(6) 



Equations (4) through (6) are translated into four pole parameter 
notation by the method of Section 2 .2 : 



n 

ro^ 




(7) 


= F 2 


station 2 blocked 


(8) 



( 
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l = z _ _ 1 

<V f 



The four pole parameters therefore become : 



(9) 



e ll 


= 1 


(10) 


£ 12 


= Z 


(11) 


e 21 


= 0 


(12) 


e 22 


= 1 


(13) 


The 


Water 





As mentioned in Section 2.4.5 the water, considered only as a 
resistance to the propagation of acoustic waves, has a specific resistance 
R* dependent on the water density and propagation velocity. This can be 
written, (2): 



R 



R* = j- 



( 1 ) 



where 



R* = Specific resistance 
R = Resistance 

A = Effective radiating area 
f = density 

c = velocity of sound in water 

Treating the water as a dashpot-type resistance (velocity 
proportional), the four pole parameters were obtained in Appendix A-2.3* 

= 1 ( 2 ) 



'll 



e l2 = ° 



'21 R 



€ 22 = 1 



(3) 

(4) 



(5) 
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R was evaluated by two methods . In one an effective area A was 
assumed constant with frequency. 

R = ^cA (6) 

The hull was considered also as a vibrating diaphragm. The radius 
r of the effective vibrating area is given by: 



d 

r = W 

where cu is the radiating frequency. 
Therefore : 



( 7 ) 



R = 




( 8 ) 



Substituting (6) or (8) into (4) gave the desired four pole para- 



meters . 

A— 5 System Properties and Systematic Changes 
A-5.1 General Description 

Four variations of the general system of Figure A-3 were analyzed 
in this study. These variations were made in the relative sizes of the 
machine and foundation structures to minimize any scale effect in the 
trends of vibration transmission. Within each of these systems, element 
properties were varied. In addition, isolation elements were inserted at 
various locations in the systems. From these variations and insertions, 
trends to provide a guide for preferred design practice in the reduction 
of vibration transmission were established. 



A-5.2 System 1 

With reference to Figure A-3 the properties of System 1 were: 
1 . Machine 

Weight of beam, W = 50*0 lbf 
Mass of beam, m = *1295 lbm 
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FIGURE A-3 
GENERAL SYSTEM 




KEY: 



1: Machiimc 

%• Moumtihg 
3: Foundation 

4 • Column 
5 * Hull 
6 : Wattr 



(6o) 



2 . 



Length of beam, L = 60.0 in 
Thickness of beam, t = .75 i n 
Mounting 

Rated natural frequency of the mount at rated load. 



U lo 

f 0 = 27= l6 '° cps 



Maximum transmission at f^, = 7*0 

Frequency of the first high frequency transmission peak. 



3- 



U t 

f . = = 1000 cps 

t 2rt * 

Rated load of mount 

Foundation 

Weight of beam, W„ 



Mass of beam, 



“b 



, W = 50.0 lbf 



= 30.0 lbf 

.078 lbm 



Length of beam, L = 36. 0 in 
Thickness of beam, t = -75 in 
4 . Column 

Circular Cross-section 
Length of column, L = 6.0 in 
Mass of column, m c = .039 lbm 

Weight of column, W c = 15-0 lbf 



5 . Hull 

The nine parameters are : 

R - 50 lbf-sec/in 
s 1 

R^ = .01 lbf-sec/in 
Rp = 5000 lbf-sec/in 
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X g = 500 rbf-sec/in 

W =10.0 lbf 
s 

g 

K g = 10 lbf / in 

W = 15 lbf 
P 

g 

K = 2 x 10 lbf/ in 
P 

6 . Water 

R = 37,600 lbf -sec/in (constant with frequency) 

A-5’3 System 2 

With reference to Figure A-3 the properties of System 2 were 
Identical to System 1 except: 

Foundation 
W n = 100 lbf 

-D 

nn, = .259 lt> m 

J3 

L = 60.0 in 
t = 1.5 in 
A~5-^ System 3 

With reference to Figure A-3 the properties of System 3 were 
Identical to System 1 except : 

Foundation 
W n = 200 lbf 

-D 

nig = .518 lbm 

L = 60.0 in 
t = 3 *00 in 
A-5.5 System 4 

With reference to Figure A-3 the properties of System k were 
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Identical to System 1 except : 



Machine 
W D = 500 lbf 

nig = 1.295 lhm 

L = 60.0 in 

t = 3.75 in 

Foundation 

= 300 lbf 
Jd 

m_, = .78 lbm 

L = 60.0 in 
t = 2.25 in 

A- 5.6 Element Insertions 
A- 5 . 6 .I Thiokol R. D. Mounting 

In the systematic study the Thiokol R. D. mounting was either 
substituted for the Portsmouth mounting or placed in parallel with it. The 
properties of the Thiokol R. D. mounting were: 

2 

Mounting cross sectional area, S = 9*^3 in 
Mounting length, L = 4.0 in 
A- 5 . 6 . 2 Foundation Damping 

Viscous foundation damping which was constant with frequency was 
first considered. The damper was inserted between the foundation and the 
column. The property of the damper was: 

System'.l : 

Damping constant, R = 390 ^ ec 



System 2 : 



R = 390 



lbf-sep. 

in 
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System 3s 



R = I56O 



System 4: 

R = 1760 



lbf-sec 

in 



lbf-sec 

in 



Damping at the same position with Thiokol R. D. as the damping 
material was next considered. The dimensions of the Thiokol R. D. were: 



2 

Cross sectional area, S =20.0 in 
Thickness, L = 1.0 in 

The thickness of the damping material was such that its compliance 
property could be neglected. 

A-5.7 Systematic Changes 

The systematic changes in the structural foundation path are best 
followed by an examination of the individual cases. 

Case 1 . System 1 used. This was the basic case for System 1 with 
which the cases with property changes could be compared. 

Case 2 ■ System 1 used. Mount damping was increased by decreasing 

Q 0 to 2.5. 

Case 3 « System 1 used. Mount stiffness was increased by setting 
= 201.062 and U = 12,566-384. 

Case 4 . System 1 used. Mount stiffness was decreased by setting 

U = 50.266 andU = 3,141.596. 

U b 

Case 5 » System 1 used. The Thiokol R. D. mounting replaced the 
Portsmouth mounting. 

Case 6. System 1 used. The hull damping was increased by setting 
R = 100, R. = .02, and R. = 2000. 

Case 7* System 1 used. The mount damping was increased by setting 
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Qq = 2.5, and. the mount stiffness was decreased by settinUJ^ = 50.2i56 
and u^. = 3,l4l.596. This case was thus a combination of cases 2 and 

4. 

Case 8 . System 1 used. The Thiokol R. D. mounting was used, and 
the hull damping was increased by setting R = 100, R = .02, and 

S X 

R^ = 1000. This case was thus a combination of cases 5 and 6. 

Case 9 - System 1 used. Mount damping was increased, mount stiffness 
decreased, and hull damping increased. This was accomplished by 
setting Q = 2.5, UJ = 50.266, U = 3.l4l.596, R = 1000, R = 1.0, 

U U o S X 

and Rg = 20,000. 

Case 10 . System 1 used. Viscous foundation damping was inserted 
in the system with R = 390. 

Case 11 . System 1 used. The Thiokol R. D. and the Portsmouth 
mounting were placed in parallel. 

Case 12 . System 1 used. The dissipative resistance of the water 

was varied with frequency. See equation (8) in Appendix A-4-5* 

Case 13 . System 1 used. The mount stiffness was decreased by u 

setting kJ = 50-266 and [J = 3,l4l-596, and viscous foundation 
0 t 

damping was inserted in the system with R = 390* This case was 
thus a combination of cases 4 and 10. 

Case l4 . System 1 used. Foundation damping with Thiokol R. D. 
was inserted in the system. 

Case 15 * System 1 used. The Thiokol R. D. and the Portsmouth 

mounting were placed in parallel. The stiffness of the Portsmouth 

mounting was decreased by setting lx) = 50.266 and U = 3>i4l.596* 

This case was thus a combination of cases 4 and 11. 

Case 16 . System 1 used. Hull damping was increased by setting 

R = 100,000, R_ = 100, and R = 2 x 10 6 . 
s * * 1 9 2 
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Case 17 . System 1 used. Foundation damping with. Thiokol R. D. 
was inserted, and the Thiokol R. D. mounting replaced the 
Portsmouth mounting. This case was thus a combination of cases 
5 and l4. 

Case l8 . System 1 used. Foundation damping with Thiokol R. D. 
was inserted, and mount stiffness was decreased by setting 
Uq = 50.266 and = 3,141.596. This case was thus a combination 
of cases 4 and l4. 

Case 19 . System 1 used. Foundation damping with Thiokol R. D. 
was inserted, mount stiffness was decreased by setting = 50.2 66 
and = 3>l4l.596, mount damping was increased by setting 
Qq = 2 . 5 . This case was thus a combination of cases 19 and 2. 

Case 20 . System 2 used. This was the basic case for System 2 
with which the cases with property changes could be compared. 

Case 21 . System 2 used. Mount stiffness was decreased by setting 
u Q = 50.206 and w = 3,141.596. 

Case 22 . System 2 used* The Thiokol R. D. mounting replaced 
the Portsmouth mounting. 

Case 23 * System 2 used. Mount damping was increased, mount 
stiffness decreased, and hull damping increased. This was 
accomplished by setting Q = 2.5* = 5 ° - 266 , \jJ- 3,1^1-596, 

U O X 

R = 1000, R = 1.0, and R = 20,000. 

S -L c. 

Case 24 . System 2 used. Viscous foundation damping was inserted 
in the system with R = 390. 

Case 25 > System 2 used. Foundation damping with Thiokol R. D. 
was inserted, and the Thiokol R. D. mounting replaced the Portsmouth 
mount ing . 

Case 26. System 2 used. Foundation damping with Thiokol R. D. 
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was inserted, and mount stiffness was decreased by setting 
= 50.266 and u^_ = 3,141.596. 

Case 30. System 3 used. This was the basic case for System 3 
with which the cases with property changes could be compared. 

Case 31. System 3 used. Mount stiffness was decreased by setting 

u» 0 - 50.266 and U t = 3,141.596. 

Case 32 . System 3 used. The Thiokol R. D. mounting replaced the 
Portsmouth mounting . 

Case 33 » System 3 used. Mount damping was increased, mount 

stiffness decreased, and hull damping increased. This was 

accomplished by setting Q = 2.5, U = 50.266, U) = 3, 141.596, 

oo t 

R = 1000, R.. = 1.0, and R. = 20,000. 

Case 34 . System 3 used. Viscous foundation damping was inserted 
in the system with R = 1560. 

Case 33 . System 3 used. Foundation damping with Thiokol R. D. 

was inserted, and the Thiokol R. D. mounting replaced the Portsmouth 

mounting. 

Case 36. System 3 used. Foundation damping with Thiokol R. D. 
was inserted, and mount stiffness was decreased by setting 
U) Q = 50.266 and = 3,141.596. 

Case 50* System 4 used. This was the basic case for System 4 
with which the cases with property changes could be compardfl. 

Case 51» System 4 used. Mount stiffness was decreased by setting 
U Q = 50.266 and (J = 3,141-596. 

Case 52 . System 4 used. The Thiokol R. D. mounting replaced the 
Portsmouth mounting. 

Case 53. System 4 used. Mount damping was increased, mount stiffness 



decreased, and hull damping increased. This was accomplished by 
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setting Q q = 2.5, W Q = 50.266, = 3,141-596, R g = 1000, R ± = 1.0, 

and = 20,000. 

Case 54. System 1 used. Viscous foundation damping was inserted 
in the system with R = I76O. 

Case 55* System 4 used. Foundation damping with Thiokol R. ti. 

i ' 

was inserted, and the Thiokol R. D. mounting replaced the Portsmouth 
mounting. 

Case 56. System 4 used. Foundation damping with Thiokol R. D. 
was inserted, and mount stiffness was decreased by setting 

lo = 50.266 and u = 3,141.596. 

U u 

Case 5$. System 4 used. Foundation damping with Thiokol R. D. 
was inserted, mount stiffness was decreased by setting = 50.266 
and = 3,141.596, and mount damping was increased by setting 

Q 0 -l 5 . 

A~6 Sample Calculations 

As noted in Section 2.5*2, the computational effort involved in 
the thesis was so great as to render it impossible without the use of a 
digital computer. However, the computational effort was largely needed 
in iterative calculations. To illustrate the complexity and length of 
the calculations involved in solving this problem, a sample calculation 
for one value of transmission is presented here. 

For example. Case 1 properties are used with the driving frequency 
at 10.3 cps or 65 rad/sec. 

The unsimplified four pole matrix equation for Case 1 follows: 
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F 1 




1 

H 

X H 

roft-< 

1 




1 


0 




M 

1 + mr 


J5J H 

1 


VI 




b 1 




i_ 


1 




M 


1 


I 








Z M 






c 





Machine Mounting Foundation 



1 Z H 




1 0 






1 


0 1 




=- 1 






R 



Hull Water 

i 

The following are the properties for this case: 

1 . Machine 



W_ = 50 lbf 

.D 

nig - .1295 lhm 

L = 60.0 in 
t = .75 in 

2 . Mounting 

^ = l6 CPS 

Q 0 = 7.0 
f - 1000 cps 

W = 50.0 

3. Foundation 
W-, = 30.0 lbf 

iJ 

m = .078 lbm 

.D 

L = 36.0 in 



t = .75 in 
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4 . Column 

L ^ 6.0 in 

m = -039 lbm 
c 

W =15.0 lbf 
c 

5- Hull 

The nine parameters are 
R g = 50 lbf-sec/in 

R^ - .01 lbf-sec/in 

R^ = 5000 lbf-sec/in 

X g = 500 lbf-sec/in 

W =10.0 lbf 



6 . Water 

R = 37,600 lbf-sec/in 

These properties are also listed in Appendix A-5-2. 

The individual four pole parameters must be evaluated before 
combination of the matrices. 

Machine : 

From Appendix A-4.1: 



s 



c 

K = 10 lbf / in 



s 



W = 15 lbf 
P 

c 

K = 2 x 10 lbf/in 
P ' 




n 



where : 
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w =71 /"§ 

n n m. 



t <§> - 



B 



^ L 
, A'n 

cosh, - — + cos g 



At 

2 

sinh/S^L + sin^^L 



cosh - COSyg^L 



...A 1 



sinh 



+ sin 



/) and Q) (tt) are listed in Reference (l6). uJ is calculated, 
n teF 2 v ' n 

= 3620 (\> ± = -1.21565 ^6 = 0 

(j) = ~1.4l4l2 

- 4.1 diooiR <t 



^ = 273 



^ - 751 



07 = 508O 



<t> 2 “ 0 



'-o = 1475 cu Q = 6780 



= 244o 



o (J 



to = 3630 to = 



(|) 3 = +1.42238 = 0 

= 0 (j> 9 = +1.41412 

d> c = -1.41386 4^0 = 0 



10 rr 10 $) 

For L0 = 65 , using the above and listed properties 

10 l 2/ L\ 

Mf Ls _ J065I S J 

2 " .1295 ; w 2 . te25 ’ (65K-1295) 



n 



M^) ^ “ -1074 j 



7T b 9>3 ° J 

M(|) 

Mounting : 

From Appendix A-4.2: 



w § 



M 



K 



1 + jQ, 



'0 



^ + 1 



^a>§ fi + — 

'* o l ti 



K = 






(71) 



For 6 5 : 



K = 



( *1295 )( loo. 53 ) 2 [l + 

" 7 Ji( 100 -537 

1 6283.2 



k te 1285 



l_ _ (6g)(7) r l + 7.1 1 

Z M 1285 l ^9 + 1 J 



=- = .00708 + .o496j 

Z M 



Foundation : 

From Appendix A-4 . 3 : 

c*=> j. /L x 



M(L) 



.i- <fc><g> 4> n < L > ,1 

m ^ P P m l i 



HU n , 1 2 , .2 

B n = 1 W -uJ 
n 



m B u^ 



The values of U) and fl) (L) were calculated as for the machine. The values 
n T n 



of(J)(^-) are the same as for the machine. 



1 

4> 



758 


1! 


10100 




= +2.0 ( 


^ - *2.0 


2085 


LUrj = 


l4o8o 




= -2.0 ( j 


l* 7 = +2 -° 


4090 


U Q = 


18780 


'fe 


= +2.0 (| 


ffe - - 2 ' 0 


6770 


II 


24200 




= -2.0 (| 


1^9 ‘ +2 -° 


10100 


0 

11 


= 30100 


4>5 


= +2.0 (j 


O 

II 

ro 

b 



5 

For U) = 65 ’ 
M(L) = 



10 



&< L > 



■1(63) \ 

^ *°^ n = 1 6J 2 - 4225 
n 



( 65 X. 078 ) 



m(l) = - .2008 j 
1 



m(lJ 



= 4.98 j 



(72) 



<*> p — 

M _ \ g( .- j-) . . i_ 

C m c n = 1 -u J 2 ) 

' n ' 

From Reference (l): 

-IV 

W = 3.33 X 10 5 
UJ 2 =9.99 X 10 5 

W = 16.65 x 10 5 

= 23.31 x 10 5 

w = 29.97 x 10 5 
0 

For id - 65 : 

5 n^l 

M = 3(65) \ 2 (-l) 2 J 

c - F039T fn (u | . te25) ^ 039,(65) 

M c = - *393 j 



Hull: 



From Appendix A-4.4 
Z = R + JX 



R = 



R_ Co 2 jtR^ R 

R + -3- + — i. + ^ 
s 2 jt 00 



R = 



1 + R 2 (^ — - ^-) 
pW _K" ' 
P P 

(101X65) . 2 jt( 1000 ) 5000 

7 2n Z3 



1 + (5TOO)2( firii5T - 



65 



2 x 10 



R = 146-704 



(73) 



OL) 

UiW K pW K ’ 
s s p p 



X = X + — + 

s g W 



1 + r 2 (S — _ JiLl. 
p'«W -K*- 
P P 



s2 ,386 



65 



x = 500 + igKiol . 10_ + 



6 w (feT - — ^ 



4- 



2 X 10 



Water : 



W 1 + ( 50°° )2( f|| 7rT5 T-^-^5-) 

X = - 14,898.3 
z = 146.7 - 14,898.3 j 

From Appendix A-4 .5 
R = 37,600 

^ = .0000266 
K 

Substituting the computed four pole parameters in equation (l) 
in polar form: 



V n 



1 9.30^/j 




1 


0 




2.96/0 


0 1 




.05/l.428 1 




-.393 


IL 


l 14,902/JS- 




1 




0 




_F 6~ 


0 .u 1 






2.66 x 10"^ / 0 


1 




V 6 






($d) 



This equation will reduce to the form: 



e ll € 12 



€ € 

21 22 



F n 




1 




V. 




1 




_ _ 





(3) 



With reference to Figure V: 



V 6 - 0 



W 



From the four poles of the element representing the water: 



( 7 ^) 



Substituting equations (4) and (5) into equation (3) : 




In section 2.5*1 the transmission T was defined as: 

T = 20 log 10 i 
Therefore : 

T - 20 lo 8 10 7 . - 

Prom obtaining equation (3)* 

€ ll - 5*51 
At 65 : 

T = - 14.85 

A-7 The IBM 709 Digital Computer Program 

j 

A- 7 .I Computer Logic 

The logical procedure comprising the program written to solve this 
problem can best be illustrated by the computer flow diagram in Figure A-4. 
A-7 *2 Scope 

The cathode ray tube oscilloscope is a convenient computer feature 
which produces a graphical representation of the problem answers. The 
computer routine used to program this feature is not included in the 
program listing. Writtenhin FAP, this routine projected a three cycle 
semi-i,og grid on the scope, labeled the grid, and plotted the answers . 
Photographs, 35 mm* in size, were taken of the scope presentations. These 
photographs appear as results in Appendix B. 



A- 7 . 3 Description of Subprograms 
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Figure A-k 



Computer Logical Flow Diagram 
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1. MAIN This serves as the control program calling all subprograms 
when needed in the solution. 

2. RDIN This subprogram reads in the data including the system 
properties. It also calls the routine to place the semi -log grid 
on the scope. 

3- MATH This program computes the four poles for the first element 
in the system. In this case the four pole parameters for the 
machine representation were obtained. 

4 . MAT 2 This serves the same purpose as MAT 1 for the second 
element or mounting in this case . 

5. MAT 3 This serves the same purpose as MAT 1 for the third 
element or foundation in this case. 

6. MAT 4 This serves the same purpose as MAT 1 for the fourth 
element or hull in this case. 

7- MAT 5 This serves the same purpose as MAT 1 for the fifth 
element or water in this case. 

8. TRANS This program takes all the four pole parameters computed 
by the various MAT routines and stores them for use in the matrix 
multiplication. It sets up the matrices prior to multiplication. 

9* MAMUL This program reduces the complex matrices to one 
resultant matrix. 

10. ANSW This program selects the matrix elements to obtain the 
answer from the final matrix of MAMUL. It also provides the routine 
for plotting the points on the scope. 

11. CHOMEG This program tests the value of frequency and trans- 
mission to insure that the scope points are numferous enough to give 
an adequate plot. It also changes u) before returning to the start 



of the program. 
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12. PRINO This program prints the values of frequency and 
transmission to conclude the program. 

A-7.4 Program Listing 

The computer program listing is included to facilitate its possible 
future use in design investigations . 
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C MAIN RRUbRAM 

DIMENSION A <23 ) *B (2 3 ) »C < 20 ) »D( 20 ) »h <400 ) *b(400) »bl 1 400 ) »V ( 200 ) . VB 
1(20) » AA ( 2 00 ) »bb ( 200 ) .CC < 200 ) »DD < 100 ) *EL ( 100 ) 

CALL KUlN ( A.B.C.o.AA.Db.CC.BB.LE) 

G( 1 ) = A ( V ) 

DO 210 L = 1 * 400 
C MATRIX CALL SEUUENCE 

CALL MAT 1 (b.A.B.C.D.L.AA.BB.CC.DD.EE) 

CALL MA r 2 (b.A.B.C.D.L.AA.BB.CC* DD * tt ) 

CALL MAT 3 ( 0 , A *B . C * D • L * AA » BB . CC »DL> . t E) 

CALL MA I 4 (b.A.B.C.D.L.AA.BB.CC.UU.EE) 

CALL MA I & ( b » A » B * C * O . L * AA * Bb * CC »DD » E t ) 

C TRANSLATION StUUtNCt 

CALL TRANS (V.EE) 

C MATRIX MOLT IRL1CAT ION 

CALL MAMOL ( 1 .J.K.V.A) 

C ANSWtRS 

CALL ANSW (A.bl.b.E.L.V) 

C CALL FRtOOtNCY CHANbt 

CALL CHOMEb (E.L.G.A) 

I F ( G ( L ) “A ( 8 ) ) 210.240.240 
210 C0NT1N0E 

240 CALL FR1NU (bl.G.E.L) 

CALL tXIT 
END 

SUBROUTINE TELE 
CALL FRAME 
CALL GRID 
CALL GRID 
CALL GRID 
CALL GRID 
CALL GRID 
CALL GRID 
CALL TITLE 
CALL VERT 
CALL HOR 
RETURN 
END 

SUBRUUTINE RD1N (A.O.C.U.AA.Bd.CC.DO.Le ) 

DIMENSION A<2s)*B<2:>)»C<20|.0<20).AA(200)»bO<200) *CC<200) .UD< lOU) . 
1 EE < 1001 

READ 213»(IA(II.1=1» 10 )»(B(I)»1=1. 20 ).(C(1I. 1=1. 20). (D(l). 1=1. 20). 
1100(11.1=1.10)) 

213 FORMAT (3E13./) 

CALL TELE 

RETURN 

END 

SUBROUTINE MAT 1 (b.A.B.C.O.L.AA.oB.CC.DO.EE) 

C Be AM MATRIX 

DIMENSION b l 400 ) »A l 2D ) *t> < 2D ) .UtO) .0(20) »AAl200) »bB l 200 ) »CC < 2UU ) . 
1DD( 100 ) »EE ( 100 ) »VB( 20 ) 

V43=1.0/(B(2)*b(L) I 
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00 206 1 = 1. 10 

IF (D ( I ) **2~0( L ) **2 ) 201*202*201 

201 VB ( I ) =0 ( L ) *C I I )**2/<D< 1 )**2-U(L)*>»2) 

GO TO 2 01? 

202 IF (G(U*C( 1 )»*2-0.0) 204*203*20A 

203 VB ( I )=0.0 
GO TO 206 

204 VB ( I ) =A ( / ) 

206 CONTINUE 

V7b=VB(l)+VB(2)+Vb(3)+VB(4)+VB(S)+Vb(6)+Vd( 7 ) +VB I 6 ) +VB ( V ) +VB ( 10) 
v/y=v7a/b i <: ) 

V229=V79~V43 

IFIV229-0.0) tub *206 *20 I 

206 V329=A( 6 ) 

GO TO 209 

207 V329= A ( 2 ) 

GO TO 209 

208 V329=0.0 
V328 = A I 1 ) 

GO TO 212 

209 V328=1.0/AB6F(V229) 

212 EE ( 1 ) = V 328 

EE ( 2 ) = V329 

RETURN 

ENO 

C THIS IS MOUNT 

SUBROUTINE MAT2 < G . A *B . C . 0 *L . AA * BB * CC *00 » EE ) 

DIMENSION 0(400) *A(26)»B(26)*C(20)*D(20)»AA(200) *BB(200) *CC(200) ♦ 
100 ( 100) *EE ( 100) 

VA = (B(4 )**2*B( 6 ) )*< (B(b )+G(L) *A( 1 ) ) / (B(6)+B(4)*A( 1) ) ) 

RE = (G(L)*B(3) )/(VA*(B(3)**2+l#0) ) 

VE = RE*B(3) 

E E ( 3 ) = SOKTF ( VE**2+RE**2 ) 

EE ( 4 ) = A T ANF ( VE / RE ) 

RETURN 

END 

SUBROUTINE MAT2(G.A.u.<-.D.L*AA*dd.CC.DD*EE) 

DIMENSION 0(400) *A(26)*bl26)*C(2U)*D(20).AA<200) *BB(200) *CC(200) ♦ 
100 ( 100) * EE l 100 ) 

GEE = 72.0*EXPF ( 0 .4* LObF ( 0 ( L ) ) ) 

RAR - 0.1*EXPF(0.326*LUoF (O(L) ) ) 

V43 = RAR** 2 
I F ( V43-1 • 0 ) 410*411.410 

410 RE = ( B( 19 ) /B( 20) )* (oEE/G(U )*( 2.0*RAR/ ( 1.0-V43 ) ) 

VA = (B(19)/B(20) )*(6EE/G(L) ) 

EE ( 3 ) = l.O/SORTF (RE**2+VA**2 ) 

EE(a) = ATANFIVA/RE) 

GO TO 412 

411 EE ( 3 ) = 0.0 
EE ( 4 ) = 0.0 

412 RETURN 
ENO 



(60) 

SUBROU I INt MA I 2 ( G *A *13 *C *D »L *AA *Ub *CC*DD*EE ) 

DIMENSION ()O.U(J ) .A< 2 3 ) *B< 2 3 ) »C( 20 ) *D< 20 ) *AA<200) *BB< 200) *CC (200) » 
1DD ( 100) * Eh ( 100) 

GEE = 72*6*EaP) (0.4*LOGMG < L) ) ) 

PAR = 0.1*EXPI < o. 325*1001- ('o(L) ) ) 

V43 = KAl<»*2 

lF(V‘+3 — 1*0) 41U»4il,«*iO 

410 RK = (B( 19 ) /D l 2 0 ) ) **■ tOtt/o(L ) )M2*0*KAR/(1»0-V43) ) 

AA = (B(19)/B120) ) M G E E / G ( L • ) 

GO TO 412 

411 HR = A ( 7 ) 

412 VA = (B< 4 ) *»2*0( ‘J ) ) * < < B <6 )+G< L ) *A ( 1 ) ) / (B (6 >+B (4 ) *A( 1 ) ) ) 

RE = (G(L)*B<3 ) )/ <VA*(B(3)**2+1.0) ) 

VE = RE*B ( 3 ) 

V47 = SUkrF(kR’"‘2+AA**2) 

V48 - ATANMAA/RR) 

V49 = SOKTr < kt 2+VE**2 ) 

V30 = ATANMVE/RE) 

TOPI = V4MV4 9 
T0P2 = V40+V30 

B0T1 = SuR I h ( ( K t f KK ) * x 2 + ( AA+VE ) **2 ) 

BO T 2 = A1ANE ( < AA » VE ) / ( RE+KK ) ) 

EE ( 3 ) = T OP 1 /BO I 1 
E E < 4 ) = TOP 2 -BO I ? 

RETURN 

END 

C THIS IS FOUNDATION 

C MATRIX THKhL 

SUBROUTINE. MA I 3 l O . A .0 * C *D *L » AA * BB ♦ CO *DD » ht ) 

DIMENSION G l 400 ) .A<23)*B<23).C(20)*D(20)*AA(200)*BB<200)*CC<200)* 
1DD( 100) »hh ( 100) * VB ( 20 ) 

V43 = 1 *0/ ( b ( /)“*>( L ) ) 

DO 303 I = 1.10 

IF<D(I+10)* tt 2-G(L)**2) 301*302.301 

301 VB ( 1 + 10 ) = (GIL )■"(_( 1 ) * C ( 1 + 10) )/(0( 1+10) *# 2 _ G ( L J **2 ) 

GO TO 303 

302 IF (G(L)*C< i )*L( l+10)-0.0) 30**, 303*304 

30 3 VB < 1 + 10) = 0.0 
GO TO 303 

304 VB ( 1+10) = A ( / ) 

303 CONTINUE 

V78 = VB ( 1 1 )+Vb ( i 2 ) + V It ( 13 )+VB ( 14 )+VB( 13 ) +VB ( 16 ) +VB ( 1 7 ) +Vb ( 18 ) + 

1 VB ( 19 )+VB(20) 

V79 = V 7 8 / 0 ( 7 I 
V22S = V79-V43 
1 F ( V229-0 « 0 ) 306*308*307 

306 EE ( 6 ) = A ( 6 ) 

GO TO 309 

307 EE ( 6 » = A ( 2 ) 

GO TO 309 

308 EE ( 6 ) = 0.0 

E E ( 5 ) = A< / ) 

GO TO 312 

309 EE ( 3 ) = 1 . 0/ABGI- ( V229 ) 
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( 81 ) 



312 RETURN 
END 

THIS 13 E Ul< EUUNDAllON WlIH DAMRLR FOLLOWING 
MATRIX fHKbt 

SUBROU1 1 NE MA I 3 (O.A.o.o.U.L.AA.HD.CL. DD » L L ) 

DIMENSION 0(400) »A( 2 3 ) • h ( 2 3 ) »(_ ( ^ U ) »D( 20 ) * A A ( 2 00 ) »IJB ( 200 ) »CC ( 200 ) * 
1DD ( 100) .ELI 100) » VO I 20 ) 

V4 3 = 1.0/ I U ( /PKjIL ) I 
DO 303 1 = 1.10 

IF(D(l+10M<*2-O(L)* ■ / ) JO 1 » 302 ♦ 301 

301 VBII+10) = (OIL t*C( I ) «L I 1 + 1 1) ) / (OI 1 + 10 ) **2-0(L ) **2 ) 

GO TO 303 

302 IF (G(L)*C( 1 )*L( I+10I-0.U) 30'*, 303. 304 

303 VB( 1 + 10 ) = 0.0 
GO TO 303 

304 VBI 1 + 10) = A ( 7 ) 

303 CONTINUE 

V78 * VBI 11 l+VB ( 12 )+VH< 13 )+VB( 14)+VB< 13 ) +VBI lb )+VB< 1 / l+VBI 18 ) + 
1VBI 19 ) +Vb I 20 ) 

V 79 = V 7B/B I t ) 

V22 9 = V/9-V43 
1FIV229-0.0) 306.308 .30 t 

306 EE I 6 ) = A ( 6 ) 

' GO TO 309 

307 EE I 6 ) = A I 2 ) 

GO TO 309 

308 E E I 6 ) = 0.0 

E E I 3 ) = A I 7 ) 

GO TO 312 

309 !EI3) = 1 .0/AB3I- I V22V I 

312 (16) = 1.0/AI10) 

ELI 14 ) = l.O+I EE I 3 ) ”LL I lb ) ) 

VA = (EE < 3 ) -"EE I lb I ) - Mm I tt I 6 ) ) 

RE = 1.0+ (EE I 3 ) »EL( lb ) ) *■ gUSE I ct (6 ) I 
EE ( 13 ) = A1ANEIVA/RL) 

RETURN 

END 

FOUNDATION Wlln uaMRlK VmKImuLl WiIm t KLUUlNl Y 
MATRIX TNKLt 

SUBROUIINE MAI3 ( o . A . B * *1 J » L » A A * am » ',.0 * uU » t_L ) 

DIMENSION Gl 400 ) »A( 23 ) »M 23 ) »C < 20) .Dl 20 ) .AA( 200 ) »I3B( 200 ) »CC ( 200 ) » 
1DD ( 100 ) .EE (100) .VB( 20 ) 

V43 = 1.0/(l)( ft 'o.j(L) ) 

DO 303 1 = 1.10 

IF(D(I+10)**2-G(L)**2) 301.302.301 

301 VB ( I + 10 ) = <G( L ) »(_ ( I ) *L ( I UO ) ) / ( u ( 1 + 10 ) * ‘2-o(L ) **2 ) 

GO TO 303 

302 IF (G(L)*C( 1 I *L< 1 + 101-0.0) 304,303.304 

303 VB ( I + 10 ) = 0.0 
GO TO 303 

304 VB ( 1 + 10) = A ( / ) 

303 CONTINUE 

V7b = Vo l 1 1 ) * vt) ( x2 ) + Vo ( L 3 ( ♦ VO ( 1 -t ) t-vo ( )j ) + V B ( 1 o ) + V ti ( 1 7 ) + VB ( 18) + 
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1VB ( 1 t )+V0 l 20 ) 

V / 9 - V/O/BI /) 

V229 = V/9-VA 3 
IFIV229-0.0) 306.30H.J07 

306 E E ( 6 ) = A ( 6 ) 

GO TO 30V 

307 EE< 6) = A ( 2 ) 

GO TO 309 

308 EE < 6 ) = 0.0 
EE ( 6 ) = A< 7) 

GO TO 312 

309 EE ( 6 ) = 1.0/ABSM V229 ) 

312 GEt = 72.6“EXRE (0.4 »LugF (G(U >) 

RAR = O.l^EXRF (0*32 6# LOBE ( G ( L ) ) ) 

V A4 = KAR«*'2 

IF (VA4-1 .0)300. 301. 300 

30G RE = B( 1 > /A ( LO ) *< GEE/G ( L ) ) * ( 2.0*RAR/ { 1.0-VAA ) ) 

Et ( 16 ) = 1.0/RE 
GO TO 302 

301 EE ( 16 ) = 0.0 

302 EE(IA) = 1.0+(EE(3)*EE( 16) ) 

VA = (EE(3)*EE(16) )*SINE(EE<6) ) 

RR = 1.0+IEh ( 3 ) *«EE ( 16 ) ) *COSF ( th(6 ) ) 

EE ( 13 ) = ATANMVA/XR) 

RETURN 

END 

C THIS IS COLUMN HULL MATRIX A 

SUBROUTINE MAI4 ( G . A . B .L . D . L . AA . BB * CC . UD . tE ) 

DIMENSION G(A00).A(23)»B(23)»C(20).D(20).AA(200) *BB(200) ,CC<200> . 
1DD ( 100 ) *EE ( 100 > » VB ( 20 ) 

VA3 = 1.0/(B(8)*G(L) ) 

DO A03 I = 1,3 

IF (DU( I )**2-G(L)**2) A01.402.A01 
AO 1 VB ( I ) = (G(L)*2»0*DD( 1+3) )/(B(8)*(DD( I ) * * 2 “G ( L ) * * 2 ) ) 

GO TO A03 
A02 VB ( I ) = A ( 7 ) 

A03 CONTINUE 

V78 = VB ( 1 ) + VB ( 3 )+VB ( 3 ) +V43 
V79 -» VB ( 2 ) + VB ( A ) 

EE ( 7 ) = SORTF (V78* # 2+V79*»2 ) 

I F ( V 7 9-0.0) AO 7,406, AO 7 
A06 EE ( 8 ) = A(6) 

A07 EE ( 8 ) = ATANF(V78/V79) 

VA5 = B ( 9 ) +(B( 1 1 ) *G ( L ) )/A(3)+A(3)*B( 12)/G(L) 

VA6 = 386.0/BI 17) *G < L ) -G ( L ) /B < 1 6 ) 

VA3 = VA3+BI 13 ) /( 1.0 + BI 13 )' 2 V46“«2) 

V4 7 = B(10)+B(1A)*G(L)/ a8b.0-B( 13)/G(L)+(B(13) **2*VA6 ) / ( 1 .0+d ( 13 ) 
1*2*VA6*#2 ) 

EE ( 9 ) = SQRTF(VA3**2+VA7»*2) 

EE ( 10 ) = AIANF ( VA//VA3 ) 

V50 = EE < 7 ) *EE ( 9 ) 

V51 = EE ( 8 ) +EE ( 10) 

V32 = V30*COSF ( VSl ) 

V53 = V30*SINF( V31 ) 
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V52 = 1.0+V52 

EE ( 1 1 ) = SQRTF(V52**2+V53**2) 

EE ( 12 ) = ATANF ( V53/V52 ) 

RETURN 

END 

MATRIX FIVE 
THIS IS WATER 

SUBROUTINE MAT 5 < G ♦ A .B . C »D .L . AA . bb ♦ CC .DD . EE ) 

DIMENSION G ( 400 ) »A(25).B(25)»C(20).D<20).AA<200) *BB(200 ) »CC<20Q) * 
1DD( 100) .EE( 100) 

EE ( 1 3 ) = 1.0/BI18) 

RETURN 

END 



SUBROUTINE MATS ( G . A .B »C . D . L .AA »BB . CC .DD ♦ EE ) 

DIMENSION GUOO) »A<25)»B<25)*C(20).D<20) »AA(200) .BBI200) #CCI200) ♦ 
1DD( 100) #EE( 100) 

EE { 13) = G(U**2/B( 1 ) 

RETURN 

END 



C MATRIX MULTIPLICATION 

SUBROUTINE MAMUL (I.J.K.V.A) 
DIMENSION V ( 200 ) .A(25) 

1 = 1 



J = 9 
K = 68 

101 K=K+2 

V41=V< I )*V< J) 

V42= V ( 1+1 )+V( J+l ) 
V43=SINF< V42)*V41 
V44=C0SF ( V42 ) *V4 1 
V45*V( I +2 ) *V ( J+4 ) 

V46 = V ( 1 +3 ) +V ( J + 5 ) 
V47=SINF(V46) *V45 
V48=C0SF< V46) *V45 
V49=V44+V48 
V50=V43+V47 

V(K)=SQRTF( V50**2+V49**2) 
1 F ( V44+V48 ) 50.25.50 

50 V ( K+l ) » ATANF ( V50/ V49 ) 



GO TO 26 

25 IF { V43+V47 ) 51» 
51 V ( K+l ) * A ( 6 ) 

GO TO 26 
27 V ( K+l ) = 

26 IFU-72) 

1100 IFIK-74) 

1101 IFIK-76) 

1102 IFIK-78) 

1103 I F < K-80 ) 

1104 IFIK-82) 

1105 I F ( K-84 ) 

1106 IFfK-86) 



27. 51 



0.0 

1099. 102. 1100 
1099. 3. 1101 

1099.4.1102 

1099.5.1103 

1099.6.1104 

1099.7.1105 

1099.8 .1106 

1099.9.1107 
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1107 1F1K-88) 

1108 IF1K-90) 

1109 IF1K-92) 

1110 I F ( K-94 ) 

1111 IF(K-96) 

1112 IFIK-98) 

1113 IFIK-1001 
1099 1 = 1+16 



1099.10.1108 

1099.11.1109 

1099.12.1110 

1099.13.1111 

1099 .14.1112 

1099.15.1113 
1099,16,1099 



J » J+16 
GO TO 101 

102 V102 = ( V< 70)*V< 72 ) )*SINF ( V( 71 )+V(73 > I 

V103 = ( V ( 70 ) *V 1 72 ( 1* COSF (V( 71 1+V( 73 ) ) 



I =■ 1-16 



J = J-14 
GO TO 101 

3 I = 1+20 

J = J+14 
GO TO 101 

4 V104 = (V(74)*V(76>>* S INF ( V ( 7 5 ) +V ( 7 7 > > 
V105 = (V(74)#V(76> )*COSF(V(7i )+V(77) > 

V ( 106 ) = V104+V102 

V ( 1 07 > = V105+V103 
I = 1-16 
J = J-16 
GO TO 101 

5 1 = 1+12 

J = J+16 
GO TO 101 

6 V102 = (V( 78>*V(80) ) *S I NF (V( 79 ) +V ( 8 1 ) ) 
V103 = (V(78)#V(80) )*COSF(V(79)+V(81) ) 

I = 1-12 

J = J— 14 
GO TO 101 

7 1 = 1+16 

J = J+14 
GO TO 101 

8 V 104 =» (V(82)*V(84) )*S1NF(V(83)+V(85) ) 
V105 » (V(82)*V(84) > *COSF <V< 83 > +V ( 8b ) ) 
V( 108 ) = V 104 +V102 

V ( 1 09 ) = V105+V103 
I = 1-20 
J = J-16 
GO TO 101 

9 I = 1+16 

J = J+1B 
GO TO 101 

10 V102 = ( V< 86) #V<86) )*S1HF ( V(tt7 )+V(89 > ) 

V 1 03 = ( V(86)*V(B8) >*CU£>F ( V(87 1+V189) ) 

1 = 1-16 

J = J-16 
GO TO 101 

11 I = 1+20 

J = J+16 
GO TO 101 

V104 = (V(90)*V(92 ) > *MNM V(9i 1+V193) ) 



12 



( 85 ) 



V 1 05 = <V(90)»V<92 I )*COSF <V(91)+V<93> ) 
v < lio) = vi04 + vi u2 
V(lll) = V109 + ViUJ 
I = I -16 
J = J-18 
GO TO 101 

13 1 = 1+12 
J = J+18 
GO TO 101 

14 V 1 02 = <V<94)*V<96> )*SINF <V(95)+V<97>> 

V 1 03 = (V ( 94 ) #V ( 96 ) ) “(.OSF ( V ( 95 ) +V (971) 

I = l- 12 

J = J-16 
GO TO 101 

15 1 = 1+16 
J = J+l 6 
GO TO 101 

16 V104 = <V(98)*V< 100) ) *6 1 NF (VI 99 ) +V < 101 ) ) 
V105 * (V(98)*V{100) )*COSF(V(99)+V<101)) 
V ( 1 1 2 ) = V104 +V102 

V ( 1 1 3 ) a V 105+V 10 3 
I = 1+85 
J = J+6 

K = ^ + 20 

17 K. = K + 2 

V43 = SQRTFIVd )**2+V< 1+1)**2) 

IF 1 V 1 1+1 ) +0*0 ) 52.26*52 

52 V44 = ATANF (V{ 1 ) / V ( 1 + 1) ) 

GO TO 29 

28 1 F ( V ( 1 ) +0 • 0 ) 53.30.53 

53 V44 = A ( 6 ) 

GO TO 29 

30 V44 = 0.0 

29 V ( K ) o V43*V< J)*SINF (V44+VI J+l ) ) 

VOC + l) » V43*V< J)*COSF< V44+VI J+l) ) 
IF(K-124)1098. 18.1115 

1115 IF1K-126) 1098.19.1116 

1116 IF (K-128 ) 1098.20.111 7 

1117 IF1K-130) 1098.21.1118 

1118 IF1K-132) 1098.22.1119 

1119 I F ( K-l 34 )1098. 23.1120 

1120 IFIK-136) 1098.24,1098 
1098 I = 1+4 

J = J+4 
GO TO 17 

18 VI 38 = V ( 122 )+V ( 124 ) 

V 1 39 = V( 123 )+V( 125 ) 

I = 1-2 

J = J-4 
GO TO 17 

19 I = 1+4 
J e J + 4 
GO TO 17 

20 V 140 = V< 126)+V< 128 ) 

V 14 1 a v< 1 2 7 ) +V ( 129) 
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I = I-S 
J = J-2 
GO TO 17 

21 1=1+4 
J = J+4 

GO TO 17 

22 V142 = V( 1 30 ) +V ( 132 ) 

V143 = V( 131 )+V(133 ) 

1 = 1-2 

J = J-4 
GO TO 17 

23 1=1+4 
J = J+4 

GO TO 17 

24 V144=V(134) +V (136) 
V145=V( 135 )+V( 137 ) 

RETURN 

END 

C THIS IS TRANSLATION 

SUBROUTINE TRANS (V.EE) 
DIMENSION V<200) »EE< 100) 
V ( 1 ) * 1.0 
V ( 2 ) = 0.0 
V ( 3 ) = EE ( 1 ) 

V ( 4 ) * EE ( 2 ) 

V ( 5 ) « 0.0 
V ( 6 ) =0.0 
V( 7) = 1.0 



V ( 8 ) 


S 


0.0 


V ( 9 ) 


3 


1.0 


V( 10) 


3 


0.0 


V( 11) 


U 


0.0 


V( 12) 


c 


0.0 


V ( 1 3 ) 


= 


EE ( 3 ) 


V( 14) 


3 


EE ( 4 ) 


V ( 15 ) 


3 


1.0 


V< 16) 


3 


0.0 


V ( 1 7 ) 


3 


1.0 


V ( 18 ) 


8 


0.0 


V< 19) 


3 


EE ( 5 ) 


V ( 20 ) 


3 


EE ( 6 ) 


V ( 2 1 ) 


8 


0.0 


V ( 22 ) 


S 


0.0 


V ( 23 ) 


a 


1.0 


V ( 24 ) 


8 


0.0 


V ( 25 ) 


3 


1.0 


V ( 26 ) 


S 


0.0 


V ( 27 ) 


3 


EE ( 9 ) 


V ( 28 ) 


3 


EE ( 10 ) 


V ( 29 ) 


a 


EE ( 7 ) 


V ( 30 ) 


3 


EE ( 8 ) 


V ( 3 1 ) 


u 


EE ( 1 1 ) ' 


V ( 32 ) 


3 


EE ( 12 ) 


VOS) 


S 


1.0 



* U VJ 
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VI 3A ) = 0.0 
V ( 35 ) » 0.0 
V< 36 ) = 0.0 
V ( 3 7 ) = EE(13I 
V ( 38 ) = 0.0 
V ( 39 ) * 1.0 
V<A0> = 0.0 
RETURN 
END 



FORTRAN 

THIS IS TRANSLATION 

THIS IS FOR FOUNDATION WITH DAMRER FOLLOWING 
SUBROUTINE TRANS 



DIMENSION VI200) .EE I 100 ) 
V ( 1 ) » l.o 
V ( 2 ) = 0.0 
V ( 3 ) » EE ( 1 ) 

VIA) « EE I 2 ) 

V I 5 ) = 0.0 
V ( 6 ) * 0.0 
V ( 7 ) = 1.0 



VI 8 ) 

V I 9 ) 
VI 10 
VI 11 
VI 12 
VI 13 
VI 1A 
VI 15 
VI 16 
VI 17 
VI 18 
VI 19 

V I 2 0 

V I 2 l 
V I 22 

V I 23 
VI 2A 
VI 25 

V I 26 

V I 27 

V I 28 

V I 29 
VI 30 

V I 3 1 

V I 32 
VI 33 
VI3A 
VI 35 
VI 36 
VI 37 
VI 38 
VI 39 
VI AO 



a 0*0 
3 1.0 
a 0.0 
= 0.0 
= 0.0 
a EEI3) 
a EE I A ) 

= 1.0 
a 0.0 
a EEIIA) 
a EE 1 15 ) 
= EE 15) 

* EE I 6 ) 
a EEI16) 
= 0.0 
= 1.0 
= 0.0 
a 1.0 
= 0.0 
a EE I 9 ) 

= EE I 10 ) 
a EE I 7 > 
a EE I 8 ) 
a EE(ll) 
a EEI12) 
= 1.0 
a 0.0 
a 0.0 
a 0.0 
a EE I 13 ) 
a 0.0 
= 1.0 
a 0.0 
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RETURN 

ENO 

SUBROUTINE ANSW ( A * G 1 ,G .E . L » V ) 

DIMENSION A(251.G1(400) »G(400) * E ( 400 ) ♦ V I 2 00 ) 

V138 = V( 122 )+V( 124) 

VI 39 =■ V< 123J+VI 125 ) 

E(L) = ( A { 5 ) / A ( 4 ) ) *LOGF ( 1 • 0/SQRTF (V138**2+V139#*2) ) 
G1(L) = G ( L > / A ( 3 > 

G2 = LOGF(GKL) )/A(4) 

MFREO = (( (G2-1. 01*920. 0) /3. 01+102.0 
MT = (( (E(L)+70. 01*870. 01/140. 01 + 152.0 
CALL BRIPICIMFREQ.MT 1 
CALL BRIP1C (MFREO. MT) 

RETURN 

END 

SUBROUTINE CHOMEG ( E .L *G. A 1 
DIMENSION E(4001 »G(400) * A ( 2 5 ) 

I F ( G ( L 1 -A ( 9 1 1251.250.251 

251 V502 = ABSF (E(L)-E(L-l 1 1 

V 1 0 1 2 = ABSF(LOGF(G(L) l-LOGF(G(L-l ) 1 1 
IF (V502-5.0 1254.253 *253 

254 IF(V502-2. 5)250. 255. 255 

255 IF(V1012-0. 028781250. 253. 250 
253 G ( L+l 1 = EXPF(LOGF(G(L) )+. 02878) 

GO TO 252 

250 G ( L + l 1 = EXPF(LOGF(G(L) 1 + . 05757) 

252 RETURN 
END 

SUBROUTINE PRINO (Gl.G.E.L) 

DIMENSION G1 (400 1 .G(400 1 .E (400 I 
PRINT 220. (G( I 1 .G1 ( I 1 »E( I 1 . I*1.L) 

220 FORMAT ( 8H OMEG = F0.1.9H FREO = F7.1.6H T = E13.6I 
CALL FRAME 
RETURN 
END 
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APPENDIX B 

PRESENTATION OF RESULTS 



The following photo plots were produced by the Cathode Ray Tube 
Recorder of the IBM 709 Digital Computer. The computed points were 
connected by solid lines by the authors to clarify the plots. 

Dhe to the space limitations of the scope presentation the 
abscissa or frequency scale is scaled in integers for three decades . 

The first decade on all of the plots extends from 10 to 100 cps, the 
second decade from 100 to 1000 cps, and the third from 1000 to 10,000 
cps . The ordinate presents transmission in decibels varying from 
~70 to +60. 
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50 lb. Machine Increased Mount Damping 

30 lb. Foundation 
Portsmouth Mounting 
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Same as Case 1, except Same as Case 1, except 

Increased Mount Stiffness Decreased Mount Stiffness 
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Decreased Mount Stiffness Replaced Portsmouth Mount 

Increased Hull Damping, Doubled 
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Decreased. Mount Stiffness 

Increased Hull Damping, 20 times original 
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Same as Case 1, except Same as Case 1, except 

Thiokol R. D. and Portsmouth Mounts Water Resistance Varied 

In Parallel With Frequency 
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Decreased Mount Stiffness 
Foundation Damping, Viscous 
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In Parallel, but Decreased Stiffness 
In Portsmouth Mount 



( 98 ) 




(99) 




(100) 




H 


to 

03 

P 0> 
p< a 


CM 


03 V» 


W 


O 

X P 


5 

O 


0) P 
CO 

•n 




O P 



<M g 
4) O 

3* 

o 

<u 

QQ 

o 3 



CO 



OQ 

a) 

V 

U 

o 

a 




CO 

ir\ 

X 

P 

0) 

83 

CM J! 

a • 

0) X>. 
•P H 

® 

SO 
CQ tf\ 






8 

•H -P 

■s § 

§ 

£ 



£ 

I 

1 

C 

£ 



( 101 ) 




Same as Case 20, except Same as Case 20, except 

Thiokol R. D. Mount Replaced Increased Mount Damping 

Portsmouth Mount Decreased Mount Stiffness 

Increased Hull Damping, 20 times 
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Foundation Dating, Viscous Thiokol R. D. Mount Replaced 

Portsmouth Mount 

Foundation Damping, Thiokol R. 
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System 3 (Appendix A-5.4) Same as Case 30, except 

50 lb. Machine Decreased Mount Stiffness 

200 lb. Foundation 
Portsmouth Mounting 



(105) 





CM 

on 

w 

o 



Thiokol R. D. Mount Replaced Increased Mount Damping 

Portsmouth Mount Decreased Mount Stiffness 

Increased Hull Damping, 20 times 
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Same as Case ju, except; 
Foundation Damping, Thiokol R. 
Decreased Mount Stiffness 
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Thiokol R. D. Mount Replaced Increased Mount uamping 

Portsmouth Mount Decreased Mount Stiffness 

Increased Hull Damping, 20 times 
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